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Dla mojej corki Natalki. 



Alice laughted. "There is no use trying" she said: "one can't believe impossible things." 
"I daresay you haven't had much practice," said the Queen. "When I was your age, 
I always did it for half-an-hour a day. Why, sometimes I've believed as many as six 

impossible things before breakfast. " 

(L. Carroll Alice's Adventures in Wonderland) 



Abstract 



We performed a search for the ''He — r] bound state via exclusive measure- 
ment of the excitation function for the dd — )■ ^Hep7r~ reaction, where the 
outgoing p — 7r~ pair originates from the conversion of the r) meson on a 
neutron inside the *^Hc nucleus. The measurements were performed at the 
Cooler Synchrotron COSY-Juelich with the WASA-at-COSY detection sys- 
tem. The internal deuteron beam of COSY was scattered on a pellet-type 
deuteron target. The data were taken during a slow acceleration of the beam 
from 2.185 GeV/c to 2.400 GeV/c crossing the kinematical threshold for the 
77 production in the dd — )■ ^Her] reaction at 2.336 GeV/c. The corresponding 
excess energy in the ^He — 77 system varied from -51.4 McV to 22 MeV. Events 
corresponding to decays of the r;-mesic ''He were selected using cuts on the 
^He momentum, p and 7r~ kinetic energies and the relative p — 7r~ angle in 
the center of mass system. The range of the applied cuts was inferred from 
simulations of decay of the '^Hc — rj bound state proceeding via excitation 
the N* resonance. The integrated luminosity in the experiment was deter- 
mined using the dd — )■ ^Hen reaction and the relative normalization of the 
points of the dd ^HepvT" excitation function was based on the quasi-elastic 
proton-proton scattering. No signal of the ^He — r) bound state was observed 
in the excitation function. An upper limit for the cross-section for the bound 
state formation and decay in the process dd {'^He — Ti)hound '^Hcp7r~, 
determined for the bound state width of 10, 20 and 30 MeV equals to 28, 32 
and 41 nb, respectively. The measured angular and momentum distributions 
of the reaction products are close to those simulated under the assumption 
of uniform phase-space distribution. 
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Chapter 1 
Introduction 



Observation of bound states of hadrons and atomic nuclei such as hyper- 
nuclei or pionic atoms opened new fields of research in nuclear physics and 
provided very fruitful results in the studies of the hadron-nucleon interac- 
tion in a many-body environment. Hypernuclei contain at least one hyperon 
in addition to nucleons. Since their discovery by Danysz and Pniewski in 
1952 [Ij, the study of their properties led to a considerable progress in un- 
derstanding the AN and interactions. In turns, pionic atoms are formed 
by a negatively charged pion trapped in the Coulomb field of the atomic 
nucleus. Observation of shifts and broadening of the energy levels in pionic 
atoms induced by the strong interaction allows for precise studies of this in- 
teraction. One of the most interesting results originating from investigations 
of deeply bound pionic atoms is the evidence for partial restoration of chiral 
symmetry in the nuclear medium [2]. 

It is also conceivable that neutral mesons such as rj, K, u, rj can form 
bound states with atomic nuclei. In this case the binding is exclusively due to 
the strong interaction and the bound state - mesic nucleus - can be considered 
as a meson moving in the mean field of the nucleons in the nucleus. Due to 
the strong attractive r^-nucleon interaction, the ?7-mesic nuclei are ones of 
the most promising candidates for such states. One expects that properties 
of ry-mesic nuclei are strongly influenced by the excitation of the iV*(1535) 
resonance which dominates the low energy rj — N interaction. The N* decays 
with roughly equal probabilities in rjN and vrA^ channel. In the second case 
it leads to the decay of the //-nucleus bound state. The discovery of the rj- 
mesic nuclei would be interesting on its own but it would be also valuable for 
investigations of the rj — N interaction and for the study of the in-medium 
properties of the N* resonance p] and of the rj meson [Ij. It could also help 
to determine the flavor singlet component of the rj wave function 0. 

The existence of ?7-mesic nuclei was postulated in 1986 by Haider and Liu 
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[0], and since then a search for such states was conducted in many experi- 
ments. However, up to now no firm experimental evidence for //-mesic nuclei 
was found. A possible reason for this can be a high background observed 
in the experiments and large predicted widths of the ?7-nucleus bound states 
ranging from 7 to 40MeV [3 El |9]. One expects that the width should be 
smaller for light nuclei, where the absorption of the t] mesons is weaker due 
to smaller number of nucleons. A very strong final state interaction (FSI) 
observed in the dd — )■ ^He?7 reaction close to kinematical threshold and inter- 
preted as possible indication of '^He — rj bound state [10] suggests, that ^He— 
system is a good candidate for experimental study of possible binding. 

Taking into account the above arguments, we proposed to perform a 
search for r^-mesic ^He by measuring the excitation function for the dd — )■ 
^HepTT" reaction in the vicinity of the t] production threshold. The outgoing 
particles in the proposed reaction correspond to the decay of ?7-mesic '^He 
proceeding via t] absorption on one of neutrons in ^He leading to excitation 
of the N* resonance which subsequently decays in the p — tt^ pair. The 
remaining three nucleons bind forming the ^He nucleus which plays a role 
of a spectator moving with relatively low momenta in the overall center of 
mass (CM) frame corresponding to the Fermi momentum distribution in 
^He. If the //-mesic ^He is produced as an intermediate state in the process 
dd — )■ ("^He — ri)hound ^Hepvr" then we expect to observe a resonance-like 
structure in the corresponding excitation function at energy below the ^B.e—7] 
threshold. 

We proposed to perform a search for the ?7-mesic ^He using the deuteron 
beam of the COSY accelerator scattered on internal deuteron target. The 
WASA-at-COSY detector was chosen for registration of the reaction prod- 
ucts. Our proposal was presented at the 34th meeting of the Program Advi- 
sory Committee (PAC) of the COSY accelerator in November 2007 [TT]. It 
received a positive opinion with the recommendation of nine days of deuteron 
beam time. The measurements were performed in June 2008. During the ex- 
perimental run the momentum of the deuteron beam was varied continuously 
within each acceleration cycle from 2.185 GeV/c to 2.400 GeV/c, crossing 
the kinematical threshold for the rj production in the dd — )■ ^Herj reaction 
at 2.336 GeV/c. This range of the beam momenta corresponds to an inter- 
val of the excess energy in the ^He — r] system from -51.4 MeV to 22 MeV. 
Unfortunately, out of nine days of the allocated beam time, due to failures 
in operation of the deuteron target and of the COSY beam we could collect 
the data only for one day. Additionally, during this period, due to the fail- 
ure of the cooling system of the WASA-at-COSY solenoid, the measurement 
was performed without the magnetic field, making impossible the momentum 
analysis of charged particles. In spite of these difficulties, the measurements 
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delivered valuable data for the search for the ^He — rj bound state. 

The present work is devoted to the analysis of these data. In Chapter 2, 
theoretical and experimental background of the search for the T^-mesic nuclei 
is presented. Chapter 3 describes the COSY accelerator and the WASA-at- 
COSY detection system. A basic concept of the experiment and simulations, 
which were performed in order to validate this concept, are presented in 
Chapter 4. Chapter 5 is devoted to the data analysis. It describes the detec- 
tor calibration, the reconstruction of dd — )■ ^Hepvr" events and the luminosity 
determination. Chapter 6 presents the final results of the data analysis in- 
cluding the excitation function for the dd — )> '^Hep7r~ process and an upper 
limit for the cross-section characterizing the production and decay of the 
77-mesic ^He. Also, analysis of experimental background, based on the recon- 
structed momentum and angular distributions of the final state particles, is 
included. Chapter 7 summarizes the thesis and provides the outlook. 



Chapter 2 

Experimental and theoretical 
background 

In the first section of this chapter we gathered basic notions and formulas 
from the scattering theory which will be used in the further parts of the 
dissertation. They concern mainly the description of bound and virtual states 
and their influence on the scattering process. For a systematic discussion 
of these topics the interested reader is referred to the cited literature. In 
the second section, the basic information about the rj meson is presented. 
The next section is devoted to the theoretical studies of the ?7-mesic nuclei. 
Also, the physical motivation of the research is presented. The last section 
contains a review of the experimental research of r;-mesic nuclei in different 
experiments. 

2.1 Bound and virtual states in the scattering 
theory 

2.1.1 Scattering matrix 

In the scattering theory the initial state of a system before collision and the 
final state after collision are connected by the scattering operator S: 

\^out)^S\^ir^)- (2.1) 

The S operator satisfies special conditions like the unitarity, the time reversal 
symmetry, the analyticity and, in the case of relativistic theory, the Lorentz 
invariance. The unitarity expresses the conservation of probability and re- 
quires that the norm of a state before and after collision has to be preserved. 
The time reversal invariance implies equality of probabilities for a direct and 
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inverse transition [12]. The analyticity properties are the consequences of 
the underlying locality of the interaction [13] . 

The initial and final state can be expanded in a basis of orthonormal 
states [T2] : 

\^out) = Y.\f){f\S\z){z\^in). (2.2) 

f,i 

The elements Sfi = {f \ S \ i) define the scattering matrix. The squared 
module l^/jp is a probability of transition from the initial state | i) to the 
final state | /). 

In the case of scattering of spinless particles the scattering matrix is 
diagonal in the basis of the angular momentum states and it depends only 
on the absolute value of the relative momentum p of the colliding particles: 

Su'ip) = Si{p)S{l~n. (2.3) 

Imposing the unitarity condition (IS'i(p)l = 1) the scattering matrix can be 
expressed in the following way: 

Slip) = e''''^^\ (2.4) 

where 6i{p) is the phase shift being a real number. The corresponding scat- 
tering amplitude equals: 

where k is the length of the wave vector equal to p/h. In the following, we 
set h = 1 and we use the relation k = p. 

The scattering amplitude expressed in terms of 6i reads: 

2i5i _ 1 1 

fl = ^ = . • (2.6) 

2tp p cot oi — ip 

2.1.2 Scattering length 

At low momenta the scattering proceeds mainly in the s-wave [l = 0) and, 
therefore, it is sufficient to take into account the //=o amplitude: 

fi=o = -. (2.7) 

p cot do — ip 

The phase shift 5o can be approximated by the effective range expansion: 

pcot5o = — + (2.8) 
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where a is the scattering length and tq is the effective range. For very small 
p: 

p cot 60 = —. (2.9) 
a 

The above equation can be considered as a definition of the scattering length, 
however, one should keep in mind, that also a different convention with 
opposite sign of the scattering length is sometimes used. With this definition, 
the scattering amplitude at low momenta can be written as: 

fi=o = = / . ■ (2.10) 

— ip — i — lap 

The scattering length is a quantity that describes the strength of the 
interaction potential in the low energy region. For repulsive potential the 
scattering length is larger than zero. For attractive interaction, positive 
scattering length corresponds to a bound state and negative one to lack of 
binding. A large negative scattering length can be connected with the so 
called virtual state which is discussed in the next subsection. 

If besides the elastic scattering, there are also inelastic channels open, 
then the scattering length has to be a complex number: 

a = aR + i ■ aj. (2-11) 

The real part of the scattering length can be interpreted as a measure of the 
elastic scattering of particles, whereas the imaginary part corresponds to the 
losses in the inelastic channels. Because of the unitarity, the imaginary part 
fulfils the condition: 

ai > 0. (2.12) 

In order to have a bound state there is a commonly quoted necessary condi- 
tion [8]: 

M > aj. (2.13) 
For a complex scattering length the fi=o amplitude has the form: 

= + . (2.14) 

The elastic cross-section is proportional to the squared scattering amplitude: 

\fi=o\' = . . , (2.15) 



From this formula it is clear that on the basis of a scattering experiment one 
can determine the absolute value of the real part of the scattering length 
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but not its sign. Therefore it is not possible to answer the question if a 
bound state exists in the system of the colliding particles. The scattering 
length can be also determined on the basis of the FSI between the produced 
particles. According to Watson and Migdal [HI [15] , the energy dependence 
of the reaction cross-section with strongly interacting particles in the final 
state can be approximated by the squared amplitude for the final state as 



given by Eq. 2.15 



2.1.3 Bound states, virtual states and resonances 

One of the most powerful techniques in the scattering theory is the analysis 
of properties of the scattering amplitude in the complex momentum or en- 
ergy plane [16]. In this approach we allow momentum to become a complex 
variable p and we treat the s-matrix elements Si{p) as analytic functions of 
complex incident momentum. The analytical properties of the s-matrix for 
given reaction depend on the properties of the potential e.g. on its asymptotic 
behavior in the limit r — > oo. For example, if the potential falls exponentially 
at large r and it is an analytic function of r in 3fJr > then the region of the 
analyticity of the s-matrix can be extended to whole complex plane except for 
finite number of singular points [16] . Alternatively, we consider the s-matrix 

2 

elements as a function of the complex energy E = However, the mapping 
from p to is a two-to-one relation. For a one-to-one relation the complex 
energy plane has to be extended to two-sheeted Riemann surface [16]. 

If the interaction is described by a real potential (no absorption chan- 
nel) then the poles can only lay either on the imaginary '^p axis or on the 
lower open half plane of the complex momenta p^. Existence of a pole of 
the scattering matrix in the complex momentum plane can be observed in 
experiment as a resonance, a bound state or a virtual state. The character 
of this physical phenomenon is determined by the position of the pole in the 



complex plane (see Fig. 2.1 ). A pole p which lays on the upper part of the '^p 



axis {^p > 0) corresponds to a bound state which is the proper eigenvalue 

2 2 

of the Hamiltonian with energy defined as = ^ = — ^ where p = ia. 
If the pole p lays on the lower half of the complex plane (Q'p < 0) then it 
can correspond to a resonance state. However, the correspondence between 
resonances and the poles in '^p < is less exact than between bound states 
and poles in '^p > 0. There can exist poles located far from the real axis that 
do not lead to observable resonance effects [TB] . The resonance state lays on 
the second Riemann sheet in the E plane and is not the proper eigenfunction 
of the Hamiltonian with the energy defined as E = ^ = Er — y. The 
resonance is represented by a non-normalizable state. 



2.1 Bound and virtual states in the scattering theory 



9 



lm{p} 



Bound state 



Re{p} 



Resonance 



Figure 2.1: Positions of poles of the scattering matrix corresponding to 
bound states and resonances. 



The simplest form of the scattering matrix which fulfils the unitarity 
condition for real p > with a pole corresponding to a bound state is: 

-p - ipo 

Si=() = : — , (2.16) 

P - Wo 

where po is real and larger than zero and p = ipo is the position of the pole on 
positive part of the imaginary axis in the complex p plane. The corresponding 
scattering amplitude equals: 

fi=o = —■ (2.17) 

-Po - W 



Comparing to Eq. 2.10 the following relation between the position of the pole 
and the scattering length can be found: 

Po = -. (2.18) 

a 

If the scattering length is much larger than the range of the potential, then 
one can derive the following relation between the binding energy and the 
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scattering length: 

-^.o.nd=|^ = ^, (2.19) 

where m is the reduced mass. This means that we can infer the binding 
energy from scattering experiments at low energies. 

A pole of the scattering matrix on the negative imaginary axis of the 
momentum plane corresponds to a virtual state. The proximity of this pole 
results in large scattering cross-sections at low energies. 

An example of an attractive interaction leading to the formation of a 
bound state is the neutron-proton interaction in the triplet state and the 
corresponding bound state is the deuteron. In turns, the neutron-proton 
scattering in the singlet state is strongly influenced by a virtual state. 

A resonance can be described by the following expression for the scatter- 
ing matrix as a function of the complex energy: 

where Efi is the resonance energy and F is the width of the resonance. The 
corresponding scattering amplitude has the form: 

and the resonance cross-section (the Breit-Wigner formula) reads: 

An (2/ + l)(r/2)^ 

In the case of several competing decay channels the spin averaged formula 
for the resonance cross-section reads as follows: 



2J+1 TT Ti-T/ 

'^'^ ' ~ {2Si + 1){2S2 + l)'f'{E- ErY + (r/2)2' ^^-^^^ 



where J is the spin of the resonance, and 5*2 are the spins of colliding 
particles, Fj and F / are the partial widths of the resonance to decay into the 
entrance and exit channels, respectively, F is the total width. 

The presence of a resonance or a bound state can be observed as a rapid 
change in the phase of S{p) and as a sharp peak in the total cross-section 
as a function of energy. As we see both, resonance and bound state are 
related. To understand the connection between resonances and bound states 
let consider the limit situation where the pole lays exactly at the boundary 
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between upper and lower half plane and let consider a potential in the form 
of XV, where A is a real coupling parameter. If we change A to obtain a more 
attractive potential then the pole moves up and becomes a bound state. If we 
make the potential less attractive then the pole will move down to the lower 
half plane C^p) and the resonance state will appear. In case of the s-wave 
limit (/ = 0) the crossing of the pole to the lower half plane corresponds to a 
"virtual state" (see Figure 2.2). Therefore, the strength of potential can have 



a capital importance in differentiating between bound states and resonances. 



Im{p} 



Bound state 



Ae{P} 



Virtual state 



Figure 2.2: Transition of the s-wave bound state into virtual state. 



The situation is more complicated if we consider also the inelastic chan- 
nels. In that case, the potential carries the imaginary part, and thus, the 
coupling parameter A is a complex number. As it is shown in [T3| the posi- 
tions of the poles in the s-matrix are shifted and the boundary limit between 
bound states and resonances is no longer the 3? axis. The poles and zeros of 



the s-matrix are "moving" as a function of the parameter A (see Figure 2.3). 
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Figure 2.3: Motion of poles of scattering matrix in the complex momentum 
plane as a function of the imaginary coupling strength A of the interaction 
potential. The idea of the picture was taken from [18]. 



2.2 Basic properties of the 77 meson 

2.2.1 Quantum numbers and decay channels 

The f] meson was discovered in 1961 by Pevsner et al. [19 jj in the Lawrence 
Radiation Laboratory. It was observed as a three pion resonance in the 
reaction 7i~^ + d ^ p+p + 7r~^ + n~ + 7r^. The experiment was performed using 
a 1.23GeV/c pion beam from the Bevatron scattered in a bubble chamber 
filled with deuterium. The trajectories of the charged pions and protons 
were measured in the chamber and the neutral pions were identified using 
the missing mass method. The t] meson was observed in the invariant mass 
spectrum of the three produced pions as a resonance with a mass of about 



550MeV/c2 (see Fig. 2.4). 



During the last 50 years since its discovery, the i] meson was intensively 
studied both experimentally and theoretically. The basic properties of this 
meson are now well established. According to the Particle Data Group [20] 
the 7] mass is equal to 547.853 ± 0.024 MeV/c^. The fuU width equals 1.3 ± 
0.07 keV and is very small compared to other mesons with higher masses. 
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Figure 2.4: Invariant mass spectrum of tliree pions as measured by Pevsner 
and collaborators [19j. Besides u meson peak at mass of around 780MeV/c^, 
a signal from rj is visible at around 550 MeV/c^. The figure was adopted from 
Ref. m. 



The values of the spin J = and an odd parity P = — 1 define rj as a 
pseudoscalar meson, rj is a neutral meson with isospin / equal to zero, even 
charge parity (C = +1) and G-parity (G = +1). 

The very small total width of the t] meson results from the fact that the 
decay into two pions is forbidden due to the parity and angular momentum 
conservation. In turns, the strong decay into three pions is forbidden by 
the G-parity conservation. Dominant are second order electromagnetic two 
gamma decay and isospin violating three pion decay (see Table 2.1). Many 
r] decay channels, which are energetically possible but are forbidden by the 
conservation of the C, P or CP symmetry are used for precise tests of these 
symmetries 
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i>lt:U.LicLl lliUU-Cb 


yi i.yu HZ u.O'rtj /o 


77 


(39.31 ± 0.20) % 




(32.57 ± 0.23) % 


Other neutral modes 


~ 0.02 % 


Charged modes 


(28.10 ± 0.34) % 




(22.74 ± 0.28) % 




(4.60 ± 0.16) % 


Other charged modes 


~ 0.76 % 



Table 2.1: Basic decay modes of 1] meson [20] . 



2.2.2 T] meson in quark model 

In the same year as the rj meson was discovered, Gell-Mann published the 
famous The Eightfold Way where he proposed classification of hadrons as- 
suming that they are built of elementary species named quarks. In his model, 
quarks appear in three fiavors: up (u), down (d) and strange (s). Barions 
consist of three quarks {qqq) and mesons are quark-anti-quark (g — q) pairs. 
Due to assumed symmetry of the strong interactions between quarks in the 
fiavor space, barions and mesons can be classified in multiplets, grouping 
particles with the same J^^ quantum numbers and with similar properties. 
According to the quark model, the rj meson can be classified as a component 
of the SU(3)-fiavour nonet of the lightest pseudoscalar mesons |22] (see Ap- 



pendix |AJ). The components of this nonet are plotted in Fig. 2.5 in terms of 
the third component of the isospin J3 and of the strangeness S. The rji and 
rjs is a fiavor singlet and fiavor octet state, respectively, with the following 
quark content: 



|?78 >= {dd + uu — 2ss) I v^, 
\rix >= {dd + uu + ss) / Vs. 



In nature, the r]i and r]s are not observed and the existing rj and rj' mesons 
are their superposition: 

\ri >= cos 6*1778 > — sin 6*1771 >, 
\ri' >= sin 6*1778 > + cos 6*1771 >, 

where 6* is the mixing angle equal to about -15° [231 El]- Because of a rela- 
tively small value of the mixing angle, the real observed rj can be treated as 
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a flavor-octet state with a small admixture of a flavor-singlet component. 



t s 



7t / 

I 

'( — 



71° 



1 K 
^ 



1 ■<\ 



-1 \ -Vi 



/ 1 



Figure 2.5: Nonet of pseudoscalar mesons ( = ). The third component 
of the isospin J3 is on the x-axis and the strangeness S is on the y-axis. 



One presumes that the flavor-singlet component r\\ can mix with pure 
gluonic states. This effect should be more pronounced in the 77' case, because 
of a greater fraction of the flavor-singlet component. However, it can be also 
important in the 77 case. According to the suggestions of Bass and Thomas 
it can signiflcantly influence properties of the ?7-meson embedded in nuclear 
matter [25] . 

2.2.3 J] nucleon interaction 

Due to the short lifetime of the r/ meson (t~ 10"^^ s) it is not feasible to 
create an 77 beam. Therefore, its interaction with nucleon or nuclei must 
be studied via the observation of flnal states of nuclear reactions including 
the ?7-nucleon (or ?7-nuclei) pair. As it will be discussed in the next section, 
the FSI between produced particles can strongly influence the production 
cross-sections and, in this way, can be used for studies of the interaction 
itself. 

In the low energy region, where the rj — N pairs are produced in the s- 
wave, the t] meson interaction with nucleons is dominated by excitation of the 
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5*11 resonance A^*(1535). The mass of this resonance is equal to 1535MeV/c^ 
and the full width is of about 150 MeV [20] ■ It is the first excited state of the 
nucleon with odd parity P = —1. As a nucleon excitation state it has the 
spin J and isospin / equal to |. Because of its large width, it influences the 
whole low energy t] — N interaction region. A^*(1535) decays predominately 
into N — 7] and N — it channels with roughly equal probability. This feature 
suggests that N — r] and N — tt pairs should be treated as strongly coupled 
systems and that the coupled channel formalism is an appropriate tool to 
describe these systems. In 1985 Bhalerao and Liu [26] performed coupled- 
channel calculations including N—rj, N—tt and A— vr channel and showed that 
s-wave 7] — N interaction is of strong and attractive nature. This result has 
raised a question whether the total interaction in a nucleus- system is strong 
enough to form a bound-state. In the bound state the t] meson undergoes 
multiple elastic scattering rjN — )■ A^* — )■ 77 — )■ A^*... until it annihilates after 
the interaction with a nucleon in the process rjN — )■ A^* — )• vrA^ (see Fig. 2.6). 




Figure 2.6: N — rj interaction presented as a series of formation and decays 
of the A^* resonance: r]N N* nN and r]N -> A^* — > r]N. The idea for 
this picture was taken from |TS] . 



2.3 Theoretical studies of ry-mesic nuclei 

2.3.1 Predictions for bounding 77 mesons in nuclei 

The existence of //-mesic nuclei was suggested by Haider and Liu in 1986 [6] . 
They investigated the ?7-nucleus interaction using the attractive 77— A^-scattering 
length ao = 0.28 -|- z0.29fm or qq = 0.27 + i0.22fm, obtained by Bhalerao 
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and Liu [26]. With this value of the scattering length r^-mesic nuclei could 
be formed for nuclei with A > 12. 

In the most of the contemporary estimations the real part of the t]N- 
scattering length is larger with its real part lying in the range from 0.5 to 
1 fm and the imaginary part equal to about 0.3 fm [27]. With larger scattering 
length, a bound state can be formed in lighter nucleus. This effect can be 
studied using the optical potential of the ?7-nucleus interaction. In the first 
order in density, this potential can be written as [28] : 

V{r) = --{l + ^)p{r)ao, (2.24) 
Ijj, rriN 

where, m^, m^v are the meson and nucleon masses, fi is the reduced meson- 
nucleus mass, ao is the r^A^-scattering length and p(r) is the nuclear density. 
Assuming fi = rriri = 547MeV/c^, rriM = 939MeV/c^, and po = 0.17fm~^ we 
obtain V=120MeV for oq = Ifm which is strongly attractive. 

Tryasuchev and Isaev performed calculations of the binding of //-meson 
in ^He and ^He nuclei using the optical potential with the nuclear density 
parametrized with the Fermi form [28] : 



Po 

1 + exp{- 



pW = ^^^^ (2-25) 



where Rc is the half-density radius, a is the thickness of the nucleus dif- 
fusion surface and po is the nucleon density of the nucleus in the center. 
They calculated the binding energy of ?7-mesons based on the exact solution 
of Schroedinger equation with the optical potential. They investigated the 
binding as a function of the real and the imaginary parts of the r/A^-scattering 
length. They concluded that for the ^He — t] system the binding is not possi- 
ble, at least for eight different values of ?7A^-scattering length taken from the 
literature. For the ^He — t] system they found that the bound state occurs in 
the case of three out of the eight considered values of the rj — A^-scattering 



length (see Fig. 2.7). These three values of the T/A^-scattering length have 



the real part larger than 0.7fm and the imaginary part is of about 0.3 fm. 



2.3.2 Motivation for a search for ?7-mesic nuclei 

With the present knowledge, the existence of ?7-mesic nuclei is not clear and 
observation of such state would be interesting on its own. In a more funda- 
mental sense, its discovery would have many important implications for the 
rj meson physics. 

The determination of the binding energy and the decay of the //-mesic 
nucleus would provide important information about the riN interaction. It 
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Re{a„} (fm) 

Figure 2.7: Curves show the boundaries of quasi-bound states in "^He — rj 
system. The darkened areas are areas of quasi-bound state formation of rj- 
meson with ^He nucleus in the complex plane of rj — A^-scattering length for 
different diffuseness a/Rc parameters: 1 - 0.25, 2 - 0.15, 3 - 0.05. Open and 
full symbols correspond to values of the ?7A^-scattering length taken from 
various works. The figure was adopted from Ref. [28]. 



would allow to verify the current values of the r^A^-scattering length which is 
only very poorly known due to the difficulties in extracting it from experi- 
mental data. 

Investigations of the //-nucleus bound states can be also important for 
studies of the properties of the iV*(1535) resonance in nuclear matter. This 
resonance dominates the r)N interaction at low energies and if its properties 
change inside nucleus, then it also influences the ?7-nucleus optical poten- 
tial [3] . The experimental investigation of T^-mesic nucleus can help to distin- 
guish between different models describing the N* structure. In chiral doublet 
models, iV*(1535) as the ffrst excited state with odd parity is treated as chiral 
partner of the nucleon. In this picture, the in-medium N* mass shift will be 
reduced due to the partial restoration of chiral symmetry. Consequently, the 
//-nucleus will have a repulsive core with the attractive part at the nuclear 
surface. On the other hand, in chiral unitary approach N* is treated as a dy- 
namically generated object in meson-baryon scattering. Here the reduction 
of the mass shift is expected to be considerably smaller and the potential 



2.4 Experimental search for ^^-mesic nuclei 



19 



stays basically attractive inside the nucleus [221 E]- 

Studies of the ?7-mesic nuclei can also help to learn about the structure of 
the 7] meson. Bass and Thomas have shown that the binding energy of the rj- 
nucleus system is sensitive to the flavor-singlet component in the rj meson [5] . 
Increasing the singlet component at the cost of the octet component results 
in a greater binding. The same happens when the gluon content contributing 
to the singlet component increases. 

The wave function of the rj meson in the bound state largely overlaps 
with the one from the nuclei. In such compact system, large medium effects 
influencing properties of mesons are expected [29j. Since the //-nucleus optical 
potential is dominated by the s-wave part, the spectroscopic studies of rj- 
mesic nuclei can provide the precise information about the s-wave potential, 
which manifest itself in the mass shift of the r] meson in the nucleus [30] . 
The estimate of the mass shift would be a very interesting result in the 
context of studies of the spontaneous chiral symmetry breaking. There are 
theoretical predictions [31j which claim that the effective restoration of the 
Uyi(l) anomalous symmetry should be observable as an i] meson mass shift 
at the flnite density in //-mesic nuclei. 



2.4 Experimental search for 77-mesic nuclei 
2.4.1 Heavy nuclei region 

Haider and Liu have shown that, with the value of the r^A^-scattering length 
determined by Bhalerao and Liu [26], the t] can be bound in nuclei with 
A > 12 [6]. Their estimations were conflrmed by the calculations of Li et 
al. [32]. Consequently, the flrst experimental efforts were concentrated on 
the heavy nuclei region. 

The flrst experiment dedicated to the search for r^-mesic nuclei was per- 
formed in 1988 at BNL [53] in nuclear reactions of the type ^y^iji^ ,p)^^^y^—rj. 
The measurements were conducted using 800MeV/c vr"*" beam and four dif- 
ferent targets: lithium, carbon, oxygen and aluminium. The momenta of the 
outgoing protons were measured with the Moby Dick spectrometer set at an 
angle of 15°. Signals from formation of //-mesic nuclei were searched in the 
energy spectra of the registered protons. In the case of formation of such nu- 
clei, e.g. in the case of the oxygen target in the reaction ^^0(7r'^,p)^^0 — rj, 
an enhancement in the proton energy spectrum at an energy close to the 
rj production threshold was expected. No such a signal was observed in 



the measured spectra (see Fig. 2.8). As suggested recently in Ref. [M], the 



kinematic conditions chosen in the BNL experiment were not optimal for 
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the production of eta mesic nuclei due to relatively large kinetic energies of 
the rj mesons with respect to the target nuclei. Therefore, further studies 
of the production of eta mesic nuclei in the (tt, A^) reaction with optimized 
kinematic conditions are planned at J- PARC |35j . 




J- 140- 
^-i 300 [■ 




220 240 260 280 300 
E (MeV) 



Figure 2.8: The proton kinetic energy spectra measured for four targets in 
the BNL search for //-mesic nuclei in reactions of the type {'K~^,p). The 
arrows indicate the eta production threshold for each target. In each case, a 
Maxwellian function was fitted to part of spectrum above the r] production 
threshold. The figure was adopted from Ref. [33] . 



Five years after the BNL experiment, a search for ?7-mesic ^^F was per- 
formed at LAMPF in pion double charge exchange (DCX) reaction ^^0(7r+, tt^) 
The 7] mesons are produced in this case in collisions of the 7r+ beam with 
neutrons inside the ^'^O nucleus: 7r~^n — )■ rjp, leading to the appearance of a 
bound ?7-^^F state. One of the possible decay channels of this state proceeds 
via absorption of the t] on one of the protons in the ^'^F nucleus, leading 
to emission of negatively charged pions: rjn — )■ 7r~p. No clear signal from 
the r] mesic ^^F nucleus was found in the measured excitation curves for the 
registered DCX process. 

In 1998 Hayano, Hirenzaki and Gillitzer [37] proposed to use the recoilless 
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((i,^He) reaction to produce t] mesic nuclei. The idea of the recoilless reaction 
is based on a choice of kinematic conditions corresponding to production 
of the T) meson at rest in the laboratory (LAB) frame so, that it can be 
easily bound by the target nucleus. The momentum of the beam is fully 
taken over by the outgoing ^He nucleus. The production of the r] mesons 
proceeds via the elementary process p{d, ^He) —t], for which the cross-section 
is relatively large even very close to threshold. Following this idea, series of 
experiments were performed at GSI, including the reactions: ^Li((i, ^He)^He— 
rj and ^^C((i, ^He)^^Be — r] |38]. Analysis of data from these measurements is 
in progress. 

Besides experimental searches for the r] mesic nuclei concluded with neg- 
ative results, there are also measurements claiming the discovery of such 
states. One of them originates from photo-production measurements on car- 
bon target performed by Sokol and Pavlyuchenko [39] . The studied reaction 
was: 

7 + ^^C ^ p{n) + i^B - T]C^C -7]) ^ n+ + n + X. (2.26) 

In this process the incoming photon produces a slow rj meson that is bound 
in the nucleus and a fast nucleon that leaves the nucleus. The decay of the 
1] mesic nucleus proceeds via the excitation of the A^*(1535) resonance and 
its subsequent decay in a vrA^ pair. In the experiment the vr+n pairs were 
registered in a two arm TOF scintillation spectrometer. The observed shift 
by 90 ± 15MeV/c^ of the invariant mass of vr+n pairs with respect to the 
A^*(1535) mass was interpreted by the authors as an effect of binding the rj 
meson in nucleus. 

Recently, the GEM-at-COSY group reported an indication of the ?7-mesic 
magnesium produced in the reaction p-|-^^Al — )■ ^He -|- ^^Mg — rj |10]. In 
their experiment, they used recoilless kinematics which means that the rj 
mesons were produced at rest with respect to the target. The momenta of the 
outgoing ^He ions were registered with the BIG KARL magnetic spectrometer 
and, additionally, 7c~p pairs originating from the t] absorption on a neutron 



were detected. Fig. 2^ shows the missing mass spectrum of the ^He ions with 
a peak around -20 MeV interpreted gnal from ^^Mg — t] production. 

In our opinion, it would be important to confirm this result with higher 
statistics. 



2.4.2 Light nuclei region 

The current estimations of the //A^-scattering length are mostly larger com- 
pared to the one used by Haider and Liu in their first predictions of r^-mesic 
nuclei. Consequently, there exist suggestions that the binding of r] mesons is 
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Figure 2.9: Missing mass spectrum converted to binding energy BE of a 
bound state ^^Mg - obtained by the GEM-at-COSY group. The figure was 
adopted from [40j. 



possible in light nuclei including ^'^He [HI US] and even deuteron [27] . In 
our opinion the light nuclei are better suited for a search of t] binding than 
the heavy ones, since due to the smaller rj absorption, the bound states are 
expected to be narrower compared to the case of the heavy nuclei. Besides, 
there are no problems with interpretation of the results connected to the 
excitation of higher nuclear levels. 

Strong indications of existence of ^He — rj and ^He — rj bound states 
originate from studies of the final state interaction in those systems. The 
early measurements of the dp — >■ ^Herj reaction close to threshold performed 
at SATURNE with the SPES-4 [B] and SPES-2 [45] spectrometers revealed a 
strong enhancement of the production cross-section due to an attractive '^He— 
r) FSI. The enhancement was interpreted by Wilkin as a possible indication 
of the ^He — rj bound state [H] . A similar effect was observed in the dd — )■ 
^He?7 cross-sections measured at SPES-4 [46] and SPES-3 [10] and it was 
interpreted as possible manifestation of r^-mesic ^He jlOj . 

Recently, high precision measurements of the dp — )■ ^He?7 reaction have 
been performed at COSY-Juelich by the COSY-11 collaboration [47] and 
independently, by the ANKE collaboration [48]. Both measurements were 
done using deuteron beam slowly accelerated in a momentum interval cover- 
ing the 7] production threshold. The total cross-sections as measured in both 
experiments are shown in Fig. 2.10 (left panel). The measurements confirm 
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a rapid increase of tlie cross-section witliin an excess energy of about 1 MeV 
above tlie threshold to a plateau of about 400 nb. The presented data points 
were parametrized with the s-wave scattering length formula [HI HTl HE] and 
a value of the ^He?7 scattering length of oq = [±(2.9 ± 0.6) + (3.2 ± OA)i] fm 
has been extracted from the fit to the COSY-11 data [47J. The performed 
analysis did not allow to answer the question if the real part of the scattering 
length is larger than the imaginary part due to a strong correlation of these 
two parameters in the performed fit. Therefore, the necessary condition for 
the existence of the bound state could not be checked. 

The COSY-11 and ANKE collaborations measured also angular distribu- 
tions of the near threshold dp — )■ ^Her/ cross-section. Close to threshold, these 
distributions can be very well described by a linear function of cos 



|.^|l + „eos9,l, (2.27) 

where a is an asymmetry parameter, which changes with energy as shown in 



Fig. 2.10 (right panel). As pointed out by Wilkin [IH], this behaviour results 
from a very strong variation of the s-wave amplitude indicating the proximity 
of a pole in the ^He — rj scattering matrix. However, information whether the 
pole lies on the bound state or the virtual part of the complex energy plane 
cannot be accessed. 

The near threshold cross-sections for the dd — )■ ^He?7 reaction are by 
about one order of magnitude smaller than the dp — )■ ^Herj cross-sections 



(see Fig. 2.11), but they show a similar energy dependence indicating for a 
strong '^He — 7] interaction in the final state. A fit of the cross-section data 
using the scattering length approximation results in a ^He — rj scattering 
length of ao = [±(3.1 ± 0.5) + ^(0.0 ± 0.5] fm This resuh can be 

converted into a pole position in the complex energy plane of \W\ ~ 4 MeV. 
Since the sign of the real part of the scattering length is not known, the pole 
position corresponds either to a virtual state or to a bound state. In the later 
case the W equals to the binding energy. 

The first direct experimental indication of a bound state of the r] me- 
son and a light nucleus was reported by the TAPS collaboration from their 
measurements of photo-production of rj mesons on ^He target [55] • Besides 
registration of 7^He — >■ ^He?7 events they measured also production of ir^p 
pairs in the process 7^He — )■ n^pX. According to their expectations, produc- 
tion of such pairs with opening angle close to 180° in the 7 — '^He center of 
momentum frame for gamma energies below the 77 threshold can indicate a 
decay of ^He — r] bound state. And indeed, the difference between excitation 
functions for two ranges of the n^p relative angle showed a structure which 



was interpreted as a possible signature of the bound state (see Fig. 2.12) 
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Figure 2.10: Left: Close-to-threshold total cross-section for the dp -^^Herj 
reaction plotted as a function of the excess energy Q. Shown are the measure- 
ments performed by ANKE collaboration (open circles) and COSY-11 
group: [HI (full dots) and [50] (triangles). The solid line represents the 
scattering length fit to the COSY-11 data [17], while the dashed line is the 
analogous fit to the data set of Ref. [H]. Right: Angular asymmetry pa- 
rameter a plotted as a function of CM momentum. Closed circles are the 
experimental data from ANKE [IH], whereas open circles represent the data 
set of COSY-11 group [1?]. The dashed and solid lines are the theoreti- 
cal parametrizations [19] explained in the text. The Figure is adapted from 
Ref. [19]. 



From a fit of a Breit-Wigner distribution to the observed structure, a bind- 
ing energy of —4.4 ± 4.2 MeV and a width of 25.6 ±6.1 MeV was deduced for 
the ?7-mesic state in ^He. As pointed out by Hanhart [SB], due to the limited 
statistics the TAPS results could also be interpreted in terms of a virtual 
^He — 7] state. Recently, the TAPS collaboration has repeated the mea- 
surements with higher statistics and after a preliminary analysis Krusche et 
al. [57] suggest that the structure observed in the n^p excitation function is 
most likely an artefact from the quasi-free 7r° production. 

A search for the '^He — rj bound state was also performed at COSY, inde- 
pendently by the COSY-TOF and COSY-11 collaboration. They measured 
the excitation function for the reaction dp — )■ pppn' close to the t] pro- 
duction threshold. One expects that pppTT~ is one of the favourable decay 
channels of the '^He — t] bound state corresponding to the r] absorption on 
the neutron inside the '^He nucleus leading to creation of the pvr" pair in the 
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Figure 2.11: Total cross-section for the dd -^^Her] reaction plotted as a 
function of the CM momentum. Shown are the measurements of Frascaria et 
al. [22] (diamonds), Willis et al. [TU] (squares), Wronska et al. [S3] (triangles) 
and Budzanowski et al. [^I] (circle). The solid line represents a fit in the 
scattering length approximation. The figure is adapted from Ref. [5^ . 
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Figure 2.12: Left and center: Excitation function of the tt" — p production for 
opening angles of 170°-180° (triangles) compared to opening angles 150°-170° 
(circles) in the 7 — ^He center of momentum system. Right: Difference of 
both distributions with a Breit-Wigner distribution plus background fitted 
to the data. The figure was adopted from Ref. [55] . 



rjn — A^*(1535) — )■ pn^ reaction . The COSY-11 collaboration observed only 
9 events which could originate from the decay of the ^He — r] bound state 
and an upper limit of 270 nb was derived for the production of such state 
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in the dp collisions [SHI EH]- The analysis of the COSY-TOF data is still in 
progress. 



Chapter 3 
Experimental setup 



The 4:71 detector facility WASA (Wide Angle Shower Apparatus) was designed 
for studies of production and decays of light mesons, especially the rare t] 
meson decays. Originally, it was installed and operated at the CELSIUS 
storage ring at the TSL in Uppsala, Sweden [60] [61], offering beams of 
protons and light ions with momenta up to 2.1 GeV/c. After CELSIUS shut 
down in 2006 the detector was moved and mounted at the COSY accelerator. 
In comparison with CELSIUS, the COSY offers higher momenta of proton 
and deuteron beams, up to 3.7 GeV/c, which allows to extend the studies 
of light mesons to a higher mass region including the (p meson. The new, 
upgraded version of the WASA detector operating at COSY is called WASA- 
at-COSY [62] [63]. 

In the first section of the present chapter basic characteristics of the 
COSY accelerator are given. The second section contains description of the 
WASA-at-COSY facility including the pellet target, the forward detector and 
the central detection system. The last two sections are devoted to the data 
acquisition system and the data analysis software. 



3.1 The COSY accelerator 

The Cooler Synchrotron C0SY[6ll [65] is an accelerator and storage ring 
equipped with phase-space cooling, operated at the Research Center Jiilich 
in Germany since 1993. A schematic view of the COSY accelerator complex 



is presented in Fig. |3.1[ The accelerator facility consists of the isochronous 
cyclotron JULIC used as an injector, the cooler synchrotron ring with a 
circumference of 184 m, as well as internal and external target stations. 
Currently, there are three internal beam experiments: ANKE, WASA-at- 
COSY and EDDA, and one detector system using extracted beam: the TOP 
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facility. 

COSY delivers beams of polarized or unpolarized protons and deuterons 
in the momentum range from 0.3 to 3.7 GeV/c. The number of stored un- 
polarized particles reaches a value of 10^^. Two beam cooling techniques, 
the electron cooling at injection energies and the stochastic cooling at higher 
energies, are applied to reduce the beam emittance and to decrease losses 
of luminosity due to heating of the beam when interacting with targets of 
the internal experiments. The typical luminosity achieved with the inter- 
nal cluster target used by the ANKE experiment is of about 10^^ cm~^ s~^ 
and with the WASA-at-COSY pellet target it is by one order of magnitude 
higher. Typical beam preparation time, including injection, accumulation 
and acceleration, is of the order of a few seconds and the beam lifetime with 
the pellet target is of the order of several minutes. One of the advantages 
of COSY is the possibility of conducting measurements during a slow accel- 
eration (ramping) of the internal beam within a given momentum interval. 
This method permits to reduce significantly a number of systematic errors 
occurring in the case when the beam is set up for each momentum separately 
and it has been successfully used in past experiments |I7] . [66] . |67] . 



3.2 The WASA-at-COSY detector 

The WASA-at-COSY detector is installed at one of the two straight sections 



of the COSY ring, right in front of the electron cooler (see Fig. 3.1). The 



WASA-at-COSY detector is depicted in the Fig. |3.2[ It consists of two main 
parts: the Forward Detector dedicated to the measurement of scattered pro- 
jectiles and target-recoils and the Central Detector optimized for registering 
of photons, electrons and pions originating from decays of mesons and excited 
baryonic states. The forward part consists of several layers of plastic scintilla- 
tors allowing for particle identification on the basis of the AE-E and AE-AE 
information and of a proportional drift chamber providing track coordinates. 
The Central Detector is composed of the electromagnetic calorimeter used 
for the energy measurement of the charged and neutral particles, the cylin- 
drical drift chamber, the superconducting solenoid providing a magnetic field 
for momentum determination of the tracks of charged particles measured in 
the drift chamber and the barrel of plastic scintillators which provides fast 
signals for the first level trigger, and together with the drift chamber and the 
calorimeter, is used for charged particle identification via AE-p and AE-E 
methods. WASA-at-COSY uses an internal target system which provides 
pellets of frozen hydrogen or deuterium. In the following subsections the 
individual components of the WASA-at-COSY detector are described. 
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Figure 3.1: Plan of the COSY accelerator complex. The beam of or, 
alternatively, D~ ions is preaccelerated in the isochronous cyclotron JULIC 
and after passing 100 m long transfer beamline, is injected into the COSY 
synchrotron storage ring via stripping injection. The operating detector sys- 
tems (ANKE, WASA-at-COSY and TOF) are highlighted in yellow. The 
places occupied by previous experiments (GEM, MOMO, HIRES, PISA and 
COSY-11) are marked as shaded. WASA-at-COSY is installed next to the 
electron cooling facility. 
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Figure 3.2: Schematic view of the WASA-at-COSY detector setup. The 
COSY beam comes from the left side. The abbreviations are explained in 
the text. 



3.2.1 Pellet target 

The WASA-at-COSY target system provides a stream of droplets (pellets) 
of frozen hydrogen or deuterium. The main parts of the system are shown 



in Fig. 3.3 



Production of the pellets starts in the pellet generator with formation of 
droplets from a jet of high purity liquid hydrogen (H2) or deuterium (D2) 
using a vibrating nozzle. The typical frequency of the nozzle vibrations is 70 
kHz. The average diameter of the droplets is around 35 /im. The droplets 
are then frozen by evaporation process in a droplet chamber. Afterwards, the 
frozen pellets enter a 7 cm long vacuum-injection capillary and obtain the 
speed of 60-80 m/s due to gas pressure difference at the ends of the capillary. 
After the coUimation process the pellets are directed through a 2 m long pipe 
into the scattering chamber and further down to the pellet beam dump. The 
basic properties of the pellet target at the interaction point are summarized 
in the Table 13. 1[ 



3.2.2 Forward Detector 

The Forward Detector (FD) provides information about charged hadrons like 
protons, deuterons or He ions scattered in the forward direction within the 
polar angle range from 2.5 to 18°. Also neutrons and charged pions can 
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Figure 3.3: Schematic view of the WASA-at-COSY peUet target system. 



Pellet diameter 


35 /im 


Pellet frequency 


8-10 kHz 


Pellet velocity 


60-80 m/s 


Pellet stream diameter at the COSY beam 


2-4 mm 


Pellet stream divergence 


0.04° 


Effective target thickness 


10^^ - 10^^ atoms cm'^ 



Table 3.1: Parameters of the pellet stream at the interaction point. 
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be measured. The FD comprises a proportional chamber of the straw tube 
type designed for tracking charged particles and a few layers of segmented 
plastic scintillation detectors used for the measurement of energy losses of 
charged particles. The particle identification in the FD is based on the AE-E 
and AE-AE information from the scintillation detectors. The signals from 
scintillators provide also the information for the first level trigger logic. The 
amount of sensitive material of 50 g cm~^ corresponds to values of 0.6 radi- 
ation lengths and 0.4 nuclear interaction lengths. The FD is placed directly 
behind the vacuum chamber, having a conical shape opening in the forward 
direction. In the forward part, the chamber contains an exit window for 
particles made of stainless steel with a thickness of 0.4 mm. The individual 
components of the FD are described in the following subsections. 



The Forward Window Counter 

The Forward Window Counter (FWC) is a thin scintillation hodoscope mounted 
directly after the vacuum chamber. It consists of two layers of 12 plastic scin- 
tillators, each 5 mm thick. The components are mounted on the paraboloidal 
stainless steel vacuum window, and they are inclined by 10° with respect to 



the plane perpendicular to the beam direction (see Fig. 3.4 (left panel)). The 
FWC is used in the first level trigger logic to reduce the background caused 
by particles scattered downstream the target. The information about the en- 
ergy loss in individual FWC layers combined with the total energy deposited 
in the Forward Detector can be used to identify charged particles and, in 
particular, to select the ^He ions. 



The Forward Proportional Chamber 

The Forward Proportional Chamber (FPC) is a tracking device placed after 
the FWC. It provides precise track coordinates (up to ~ 0.2° angular reso- 
lution [691) of charged particles passing through. It consists of 4 quadratic 
modules with a circular opening for the beam pipe at the center. Each 
module has 4 layers of 122 proportional straw tubes. For a 3-dimensional 
reconstruction of multi-track events the straw tubes in consecutive modules 
are oriented at +45°, -45°, 0° and 90° with respect to the vertical direction 



(see Fig. 3.5). The straws have 8 mm diameter and are made of thin (25 /xm) 
mylar foil coated with 0.1 fim aluminium on the inner side only. The 20 /im 
diameter anode wire made of gold plated tungsten is placed in the center of 
each straw. The FPC works with Argon-Ethane 20%-80% gas mixture at 
atmospheric pressure. 
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Figure 3.4: Front view of a one layer of Forward Window Counter containing 
12 segments of scintillation detectors (right panel) and the side view of the 
detector (left panel). The picture is taken from [68] . 



The Forward Trigger Hodoscope 

The Forward Trigger Hodoscope (FTH) is made up of three layers of 5 mm 
thick plastic scintillators. The first and the second layer consist of 24 Archi- 
median spiral shaped segments oriented clock-wise and counterclock-wise, 
respectively. The third layer has 48 scintillator segments with a torte-like 



shape (see Fig. 3.6, left part). The overlap of hit segments in the three layers 



allows to localize tracks with the FTH (see right part of Fig. 3.6). This pro- 
vides fast information about the track polar angle for the triggering system. 
The FTH supplies also the information about the track multiplicity and the 
energy losses for the trigger. 

The Forward Range Hodoscope 

The Forward Range Hodoscope (FRH), which is situated behind the FTH, 
consists of five layers with 24 plastic scintillator elements each. The first 
three layers are 11 cm thick, while the scintillator elements in last two layers 
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Figure 3.5: Forward Proportional Chamber 3-D layout. 



are 15 cm thick (see Fig. 3.7). The information from FRH combined with 



FTH and FWC is used for energy determination of charged particles and 
for particle identification by means of AE-E method. The reconstruction of 
the kinetic energy and the identification of charged particle are based on the 
pattern of deposited energy in the different detector layers. The maximum 



kinetic energy for particles stopping in the FRH is listed in Table 3.2 The 



energy resolution for protons, deuterons, and alpha particles stopped in the 
detector is approximately 3%. 



The Forward Range Intermediate Hodoscope 

The Forward Range Intermediate Hodoscope is an additional scintillator 
hodoscope, which provides two-dimensional position sensitivity. It can be 
mounted between the third and fourth layer of the FRH. This detector was 
not used during the experiment reported in this dissertation. 
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Figure 3.7: Scheme of the Forward Range Hodoscope. The specified diam- 
eter of the 1st and 4th layer is given in [mm]. 
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Particle 


Maximum stopping energy [MeV] 




200 


P 


360 


d 


450 




1000 


^He 


1100 



Table 3.2: Maximum stopping energies in the Forward Range Hodoscope. 



The Forward Veto Hodoscope 

The Forward Veto Hodoscope (FVH) is the last detector layer of the FD. 
It consists of twelve horizontal plastic scintillator bars equipped with photo- 
multipliers on both sides forming a wall (see Fig. 3.8). The information from 
the FVH is used in the first level trigger logic to select or reject particle that 
punched through the FRH. 




Figure 3.8: Scheme of the Forward Veto Hodoscope. 
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The Forward Absorber 

The Forward Absorber (FRA) is a layer of iron which can be introduced 
between the last layer of FRH and the FVH, to discriminate the fast protons 
from elastic scattering from slower protons originating from investigated re- 
actions. The thickness of the absorber can be chosen between 5 cm and 10 
cm. It has been used for example in the pp — )■ ppr] measurement. In this 
case, the fast protons from clastic scattering penetrated FRA, whereas the 
protons associated with the rj production were stopped. This absorber was 
not used during the present experiment. 

3.2.3 Central Detector 

The Central Detector (CD) surrounds the interaction point and is designed 
mainly for detection and identification of decay products of 7r° and rj mesons: 
photons, electrons and charged pions. The CD consists of the Mini Drift 
Chamber, the Superconducting Solenoid, the Plastic Scintillator Barrel and 
the Scintillation Electromagnetic Calorimeter. 

The Superconducting Solenoid 

The Superconducting Solenoid (SCS) is installed inside the calorimeter and 
encloses the MDC and the PSB detectors. It provides an axial magnetic 
field along the beam for the momentum determination of tracks of charged 
particles measured in the MDC. It shields also the detector parts against the 
high flux of low energy delta electrons produced in the interaction region. The 
SCS provides a field of l.STesla at the interaction point. In order to reduce 
the probability of electromagnetic showers caused by gamma conversion in 
the material, the wall thickness of SCS is minimized to only 0.18 radiation 
lengths. The flux of the magnetic field outside the solenoid is closed by 
means of an iron yoke. The yoke serves also as support for the calorimeter 
crystals. In normal operating mode the SCS is cooled using liquid helium to 
the temperature ~4.5 K. 

During the present experiment the cooling system of the SCS was broken 
and, therefore, no magnetic field was provided. Thus, the momentum analysis 
of charged particles registered in the MDC was not possible. 

The Mini Drift Chamber 

The Mini Drift Chamber (MDC) surrounds the beam-target interaction re- 
gion. It is used for the determination of charged particle momenta and the 
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reaction vertex. It covers scattering angles from 24° to 159°. The angle res- 
olution provided by the MDC is about 1.2°. The MDC is composed of 17 
cylindrical layers with 1738 straw tubes in total. The diameter of the straw 
tubes in the 5 inner layers is 4 mm, 6 mm in the 6 intermediate layers and 
8 mm in the 6 outer layers. The straws are made of 25 /im mylar foil coated 
with 0.1 /zm aluminium on the inner side. A 20 /im diameter anode wire made 
of gold plated tungsten is placed in the center of each straw. The straws in 
the nine inner layers are parallel to the beam axis (z-axis). The next 8 layers 
have small skew angles (6-9°) with respect to the z-axis. These layers form 
a hyperboloidal shape. The MDC is fitted inside a cylindrical cover made 
of 1 mm thick Al-Be. The straws in each layer are mounted between 5 mm 
thick Al-Be end-plates. The layers are assembled around 60 mm diameter 
beryllium beam pipe. The wall thickness of the beam pipe is 1.2 mm (see 



Fig.g. The MDC works with Argon-COa 50%-50% gas mixture. The front 
end electronics for the MDC is based on the CMP16 amplifier-discriminator 
chip originally developed for the CMS experiment at CERN ^70j|. A detailed 
description of the MDC can be found in [71] . 




Figure 3.9: Left: The assembled MDC inside the Al-Be cylinder. Right: 
Cross formed by the beam pipe (vertical element) and part of the target pipe 
(horizontal tubes). 



The Plastic Scintillator Barrel 

The Plastic Scintillator Barrel (PSB) is located inside the SCS coil and sur- 
rounds the MDC. It provides fast signals for the first level trigger and together 
with the MDC and the calorimeter it is used for charged particle identifica- 
tion by the AE-p and AE-E method, respectively. It serves also as a veto 
counter for 7 identification. In total PSB contains 146 elements of fast plas- 
tic scintillator, each 8 mm thick. PSB is composed of cylindrical part and 
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two end-caps. The cyhndrical part consists of 50 scintiUator bars arranged 



in two layers. The bars are 550 mm long and 38 mm wide (see Fig. 3.10). 
The neighbouring bars overlap by 6 mm to avoid that particles pass without 
registration. Two bars in the top and in the bottom part of the cylinder are 
split to leave space for the target tube. The end-caps contain 48 " cake-piece" 
shaped elements each. The front end-cap is flat while the rear end-cap forms 
conical surface. Both end-caps have a central hole for the beam pipe. Scin- 
tillators are glued to light guides coupled to the photomultiplier tubes. The 
photomultipliers are installed outside of the iron yoke to shield them against 
the magnetic field. The length of the light-guides is about 500 mm. 




Figure 3.10: Schematic view of the Mini Drift Chamber (brown) enclosed 
by the Plastic Scintillator Barrel (blue). On the top part, the hole for the 
target tube is visible. 



The Scintillation Electromagnetic Calorimeter 

The Scintillation Electromagnetic Calorimeter (SEC) is placed outside the 
SCS magnet. It is used to measure the energy of charged and neutral par- 
ticles in the CD. The energy resolution for 0.1 GeV photons is about of 8% 
and for stopped charged particles is about 3%. The SEC can provide also an- 
gular information with a scattering angle resolution of about 5°. The energy 
threshold for detection of photons is about 2 MeV. The basic SEC parameters 



are presented in the Table 3.3 
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Amount of sensitive material 


135 g cm ^ 


[radiation lengths] 


16 


[nuclear interaction lengths] 


0.8 


Geometric coverage (Air) 


96% 


polar angle 


20°-169° 


azimuthal angle 


0°-360° 


Maximum stopping kinetic energy 




7r±/p/d 


190/400/450 MeV 


Angular resolution 


5° 


Time resolution 




photons 


40 ns 


charged particles 


5 ns 


Relative energy resolution (FWHM) 




photons (0.1 GeV) 


8% 


stopped charged particles 


3% 



Table 3.3: Basic parameters of the Scintillator Electromagnetic Calorimeter. 



The SEC is composed of 1012 sodium-doped Csl scintillating crystals 
placed between the superconducting solenoid and the iron yoke. It covers 
the scattering angles in the range from 20° to 169° and nearly 100% of the 
azimuthal angle providing about 96% of geometrical acceptance. The crystals 
are placed in 24 layers along the beam. There are three main parts in the 
SEC: forward, central and backward. The forward part covers a scattering 
angle range from 20° to 36° and consists of 4 layers with 36 elements each. 
The central part covers scattering angle range from 36° to 150° and consists 
of 17 layers with 48 elements each. The backward part covers the region 
from 150° to 169° and consists of 3 layers. The layer closest to the beam 
pipe has 12 elements and the two outer layers consists of 24 elements each. 



The angular coverage of the SEC is presented in the Fig. 3.11 The crystals 



have the shape of a truncated pyramid. The length of the crystals varies 
from 30 cm (central) to 25 cm (forward) and 20 cm (backward part). The 
length of the crystals corresponds to the value of ~16 radiation lengths and 
~0.8 nuclear interaction lengths. The crystals are connected by plastic light 
guides with the photomultipliers placed outside of the iron yoke. More details 
about the SEC can be found in [72] . 
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ej.b(deg) 

Figure 3.11: (Top) Cross-section of Scintillation Electromagnetic Calorime- 
ter. (Bottom) The angular coverage of the SEC. The 17 layers of central part 
are marked in bright colour. The first three layers on the left (dark colour) 
belong to the backward part and the last 4 layers on the right (dark colour) 
belong to the forward part. The length of the crystals are marked on the 
y-axis. The number of elements in appropriate layer is given above the plot. 
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3.3 Data Acquisition System 

The Data Acquisition (DAQ) for the WASA-at-COSY detector is based on 
the third generation of the DAQ systems used in experiments at COSY 



The overview of the DAQ is presented in Fig. 3.12 
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Figure 3.12: Schematic overview of the data acquisition system. The figure 
is taken from [74|| . 



The DAQ is based on FPGA-controlled read out boards used for digiti- 
zation and buffering of data. It permits to conduct measurements at average 
accepted event rate of ~10000 1/s [75J with a life time of 80% of the sys- 
tem. The limiting factors are the events size and the writing speed to disk 
storage of approximately 80MB/s. The DAQ system runs in "common stop 
mode" which means that the trigger signal is coming after the data have 
been digitized. 

The analogue and digital signals from the detector front-end cards are dig- 
itized by means of Charge-to-Digital Converter (QDC) and Time-to-Digital 
Converter (TDC) read out boards. There are two types of the QDC boards: 
SlowQCD and FastQDC as well as two types of TDC boards: SlowTDC and 
FastTDC. The SlowQDC is a 16 channel Flash ADC board designed for the 
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readout of the SEC photomultipliers. It provides time stamps along with 
the charge integration. The sampling frequency is 80 MHz. The FastQDC 
is a 16 channel Flash ADC board with the higher sampling frequency of 160 
MHz. It is designed for the readout of the plastic scintillator detectors. The 
SlowTDC is a 64 channel Fl ASIC board [76]. It is designed for the readout 
of straw tubes from the MDC and the FPC detectors. Finally, the FastTDC 
is a 64 channel GPX ASIC board designed for time stamping of the plastic 
scintillator signals [77] . 

The synchronization system is used to control and synchronize the data 
flow. When the trigger electronics generates a trigger signal, the synchroniza- 
tion system issues an event number with a time stamp, which is distributed 
to all QDC and TDC boards. 

Read out boards collect and store the digitized signals and mark them 
with time stamps. They work in self-triggering mode. When the trigger 
arrives, the signals inside a given time interval are selected based on the 
trigger's time stamp. The data are then sent via high speed optical links 
to the PC computer farm and further to the event builder. The events are 
flnally written to the disk storage. 

More information about the DAQ system can be found in [73] . 

3.4 Analysis software 

For simulations of measurements with the WASA-at-COSY detector Wasa 
Monte-Carlo (WMC) package based on the GeantB software has been devel- 
oped. In turns the data analysis is performed with the RootSorter which is 
based on the Root framework. Both packages work in the Linux environment. 



Chapter 4 

Concept of experiment 



Present chapter describes the basic idea of our experimental search for the 
^He — rj bound state and results of simulations which were performed to 
validate this idea. 



4.1 Basic idea 

In our experimental studies, we use the deuteron-deuteron collisions at ener- 
gies below the r] production threshold for production of the r] — "^He bound 
state. We expect, that the decay of such state proceeds via absorption of 
the 7] meson on one of the nucleons in the ^He nucleus leading to excitation 
of the A^*(1535) resonance which subsequently decays in pion-nucleon pair. 
The remaining three nucleons play a role of spectators and they are likely to 
bind forming ^He or nucleus. This scenario is schematically presented in 



the Fig. 4.1 



According to the discussed scheme, there exist four equivalent decay chan- 
nels of the C^He — ri)hound state: 



(^He - r])bound ^Hepvr 
(^He - r])bound ^Hen7r° 
(^He - v)bound ^ ^^HpvrO 
(^He - ri)bound -> ^H?27r+ 



In our experiment we concentrated on the first one out of the listed de- 
cay modes due to the highest acceptance of the WASA-at-COSY detector 
in this case. The outgoing ^He nucleus plays the role of a spectator and, 
therefore, we expect that its momentum in the CM frame is relatively low 
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Figure 4.1: Schematic picture of the (^He — ri)bound — ^ ^Hepvr" decay. In 
the first step the t] meson is absorbed on one of the neutrons and the N* 
resonance is formed. Next, the N* decays into a p — n~ pair. The ^He plays 
the role of a spectator. 



and can be described by the Fermi momentum distribution of nucleons in the 
^He nucleus. This signature allows to suppress background from reactions 
leading to the ^Hepvr" final state but proceeding without formation of the 
intermediate (*^He — ri)bound state and, therefore, resulting on the average in 
much higher CM momenta of ^He. A kinematic variable correlated with the 
^He CM momentum is the relative angle of the outgoing p — pair. In the 
limit of ^He produced at rest in the CM frame this angle is exactly equal 
to 180° but due to the presence of the Fermi motion it is smeared by about 
30°. Quantitative predictions for the ^He CM momenta and for the p — 
relative angle are given in the next section. 

The principle of the present experiment is based on the measurement 
of the excitation function of the dd — )■ ^Hepvr" reaction for energies in the 
vicinity of the t] production threshold and on the selection of events with 
low ^He CM momenta. In the case of existence of the ^He — t] bound state 
we expect to observe a resonance-like structure in the excitation function at 
CM energies below the t] threshold. From the central energy of the observed 
structure E'-^^ one can determine the binding energy of the (^He — r/) system: 

Ebe = mne + - E^^' . (4.1) 

The width of the structure is equal to the width of the bound state. 
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4.2 Simulations of the dd — > (^He — ri)hound 
^He^^TT" process 

In order to check the feasibihty of the search for the ^He?7 bound state with 
the WASA-at-COSY detector we performed simulations of the dd — )■ ^He — 
v)bound — ^ ^Hepvr" process. These simulations were also used in preparatory 
phase of the experiment for setting adequate triggering conditions of the 
detector read out and, further on, during the data analysis for choosing 
optimal selection criteria and cuts. We assumed that the decay of ^He — 1] 
bound state proceeds via absorption of the rj meson on one of neutrons in 
the '^He nucleus leading to excitation of the A^* resonance which subsequently 
decays into the p — n~ pair. The remaining three nucleons bind forming the 
^He nucleus. In the simulation, the direction of the '^He momentum vector 
is selected assuming its isotropic distribution in the CM frame. The length 
of the momentum vector is randomized according to the probability density 
distribution of the Fermi momentum in '^He taken from [78l based on 1791: 



fip) 



ai 



■ exp "2 



(4.2) 



where the values of the parameters are: ai = 0.0001989184519 (GeV/c)'^, 
02 = 0.0028615450879 (GeV/c)^. This distribution has a maximum at 
around p = 0.18GeV/c (see Fig. 4.2). 
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Figure 4.2: Probability density distribution of the Fermi momentum in the 
^He nucleus taken from [ 78 j. 
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The four-momentum of the intermediate N* state is calculated using the 
energy-momentum conservation principle: 

PN* = Ptotal - PHe = (v^ - ^^f , -pgf )• (4.3) 

The simulation of the N* decay into a pTT~ pair is performed under as- 
sumption of its isotropic angular distribution in the A^* rest frame. The 
opening angle between the outgoing proton and pion is equal to 180° in the 
A^* reference frame, whereas in the CM frame its maximum is shifted to 
smaller values and it is smeared by about 30° due to Fermi motion of the 



nucleons inside the ^He nucleus (see Fig. 4.3). 

As a model for background processes in the present studies we consider 
the reaction dd — )■ ^Hepvr" for which we assume a uniform distribution of the 
reaction products over the available phase space. Further on, we refer to this 
process as a direct production. In contrast to the dd — )■ i^B.eri)hound ^B.epTT~ 
reaction, in the direct production the distribution of the relative proton-pion 
angle in the CM system 0^*^ covers the full angular range (see Fig. |4.3). 



The relative p — 7r~ angle is strongly correlated with the ^He mo- 

mentum p"^^ . In the case of ^He produced at rest in the CM system, 0^*^ 
equals 180°. The Fermi motion results in smearing of both the Qpl!^ and 
P%^ ■ -^o^ direct production the p^^ distribution is much wider than for 



the He — rj bound state decay (see Fig. 4.4). As it will be discussed in the 
next chapter devoted to the data analysis, a cut on the ^He momentum in 
the CM frame is used as a basic mean to suppress the background. 

For preparation of experimental trigger we performed estimation of angu- 
lar acceptance of the WASA-at-COSY for registration of the decay products 
of the '^He — rj bound state. As shown in the Fig. 4^ most of the ^He 
ions (81 %) are emitted in the forward direction and are confined in the 
acceptance of the Forward Detector. In turns, the outgoing protons and neg- 
atively charged pions are registered at much larger angles and about 71% of 



all p — Tc~ pairs is registered it the Central Detector (see Fig. 4.6). There- 
fore, in the experiment the main trigger dedicated to the registration of the 
dd — )■ C^He — rj) bound — ^ ^B.ep7i~ process was constructed under assumption 
that the ^He ions are detected in the Forward Detector and the p — n pairs 
are registered in the Central Detector. The overall geometrical acceptance is 
equal to about 60%. 



4.2 Simulations of the dd — t- ("^He — r])bound ^Hepir process 



48 
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Figure 4.3: Distribution of the p — 7r~ opening angle in the CM system 
obtained in simulation of the processes leading to the creation of the ^He?7 
bound state: dd — )■ C^Herj) hound ^Hepvr" (red line) and of the direct dd — )■ 
^Hepvr" reaction (black line). The simulation was done for momentum of the 
deuteron beam of 2.307 GeV/c. 
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Figure 4.4: Distribution of the ^He momentum in the CM system obtained 
in simulation of the processes leading to the creation of the '^B.er] bound state: 
dd — )■ {^tleri)bound — ^ ^B.epn~ (red line) and of the direct dd — )■ ^Hepvr" decay 
(black line). The simulation was done for momentum of the deuteron beam 
of 2.307 GeV/c. 
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Figure 4.5: Distribution of the ^He scattering angle seen in the LAB frame as 
simulated for the processes leading to the creation of the ^He — r] bound state: 
dd —J- i^Her]) bound — ^ ^Hepvr". The red area represents the WASA-at-COSY 
Forward Detector acceptance. 




Figure 4.6: Distribution of p (x-axis) versus tt~ (y-axis) scattering angle 
seen in the LAB frame as simulated for the processes leading to the cre- 
ation of the ^He — t] bound state: dd — )■ {'^Her]) hound ~^ ^Hepn^ . The black 
rectangle corresponds to the range of the WASA-at-COSY Central Detector 
acceptance. 



Chapter 5 
Data analysis 



The present search for the ^He — 1] bound state was performed in a dedicated 
experiment in June 2008. During the measurements, a scan of the excitation 
function for the process dd — )■ (^He — ri)bound — ^ ^Hepn^ was conducted using 
slow acceleration of the COSY deuteron beam. The integrated luminosity 
in the experiment was determined using the dd — )■ ^Hen reaction and the 
relative normalization of the data points measured as a function of the beam 
momentum was based on the quasi-elastic proton-proton scattering, which 
was registered with relatively high statistics. 

The basic settings of the COSY acceleration cycle and the trigger applied 
for the detector read out are presented in the first section of the present 
chapter. The calibration of the WASA-at-COSY detectors is presented in 
the second section. The main part of the data analysis is devoted to selection 
and reconstruction of the dd ^Repn' events with subsequent application of 
cuts reducing the counts not originating from the searched ('^He — rj) bound — ^ 
^Hepvr" decay. This part of the analysis is described in the third section. 
The luminosity determination is presented in the section four. 

5.1 Accelerator cycle and experimental trig- 
ger settings 

The presented experiment was planned for 9 days of the beam-time which 
were granted by the COSY Programme Advisory Committee, however, due 
to severe problems with the COSY beam and with the WASA-at-COSY 
deuteron target, the data was taken during one day only. The analysed 
data set consists of 66 runs (run numbers: 9163-9228) and corresponds to an 
effective measurement time of about 16.5 hours. 



Setting of the accelerator cycle is specified in Tab. 5.1 Duration of the 
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cycle was 120 s. In the first 3.431s the beam was accelerated in a routine 
way during a fast ramping of the COSY dipole magnets to the momentum 
of 2.185 GeV/c. After this, a slow ramping phase taking 107.6 s followed. 
During this phase the beam momentum was increased linearly in time to 
upper limit of 2.400 GeV/c. Events containing the detector response were 
stored together with a content of a precise clock representing the informa- 
tion about the instantaneous beam momentum. The data was taken during 
the slow acceleration phase in the momentum range from 2.192 GeV/c to 
2.400 GeV/c. During the data analysis, this range was divided into 20 inter- 
vals of equal width and the excitation function was determined on the basis 
of events selected in each interval. 



beam cycle time 


118 s 


start DAQ 


7 s 


start slow ramping 


3.431 s 


slow ramping time 


107.6 s 


lowest beam momentum 


2.185 GcV/c 


highest beam momentum 


2.400 GeV/c 


effective measurement time in the cycle 


88% 



Table 5.1: Setting of accelerator cycle. 



Several hardware triggers were used in the experiment for the read out of 
events. The main trigger dedicated to the study of the dd — )■ C^tie—rj) bound — > 
^HepTT" reaction required at least two charged tracks in the Central Detector 
for registration of the p — 7T~ pairs and at least one charged track in the 
Forward Detector for registration of the "^He ions. Additionally, a high energy 
threshold was set in the first FWC layer and the first FRH layer in order to 
suppress signals from fast protons and deuterons in the Forward Detector. 
Besides, in order to increase a selectivity for detection of ^He ions which are 
expected to stop in the first two layers of the FRH, a veto was set on the 
third FRH layer. 

For the luminosity monitoring two additional triggers were used. The 
dd ^Hen reaction was registered for determination of absolute value of 
the integrated luminosity. The measurement was based on detection of the 
outgoing ^He ions in the Forward Detector. Therefore, the applied hardware 
trigger required at least one charged particle in the Forward Detector and, 
in addition, a high energy threshold in the FWC. The prescaling factor for 
this trigger was equal to t^. For determination of integrated luminosity 



5.2 Detector calibration 



52 



as a function of the beam momentum the quasi-elastic p — p scattering was 
measured. For this, at least one charged particle in the Forward Detector 
and also at least one charged track in the Central Detector was required. 
The prescaling factor was equal to 

5.2 Detector calibration 

Registration of the reactions dd — )■ ^Hepn~ and dd '^Hen, which are the 
main processes measured in the present experiment, is based on the energy 
deposits of the ^He ions in the Forward Range Hodoscope. Therefore, a 
precise energy calibration of the FRH is of high importance for the analysis of 
theses channels and, in particular, for a clean identification of the '^He ions on 
the basis of the energy losses in the FRH layers. The calibration of the FRH 
is described in details in the first subsection. The second subsection presents 
the energy calibration procedure of the Electromagnetic Calorimeter which 
was used for identification of the protons and pions from the dd — > ^Hepvr" 
reaction on the basis of the energy losses. The position calibration of the 
Straw Tube Chambers (FPC and MDC) which were used for reconstruction of 
trajectories of charged ejectiles including ^He, protons and pions is described 
in the third subsection. The calibration for other detectors and are not 
discussed here. Their settings were taken from the previous measurements. 

5.2.1 Forward Range Hodoscope 

The calibration of the plastic scintillator detectors in the FRH comprises 
translation from the ADC channels to the energy. The calibration is per- 
formed in two steps. First, a deviation from the uniform light collection 
efficiency over the area of the scintillator detector is determined and cor- 
responding corrections to the ADC channels are calculated. In the second 
step, the conversion of the corrected ADC values to the deposited energy is 
performed. 

The non- uniformity of the light collection is caused by the fact that the 
photomultipliers reading out the scintillator segments in the FRH are at- 
tached to these segments on the outer rim of the FRH layers. Therefore, 
particles scattered at small angles interact with the segments close to the cen- 
ter of the layers, and the produced scintillation light has to pass the longest 
way to the photomultipliers. This leads to the largest light losses. In turns 
the particles scattered at the maximum angle accepted by the FRH (18°) pro- 
duce scintillations close to the photomultipliers and thus the light losses are 
relatively small. The non-uniformity of the light collection is parametrized as 
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a function of the scattering angle 6. To correct for the non-uniformity effect 
we select fast particles, such as protons from the quasi-elastic scattering, de- 
manding exactly one reconstructed charged track reaching the fifth layer of 
the FRH. These particles are expected to be close to the ionization minimum 
which means that the energy deposit per unit length does not depend on the 
scattering angle. The ADC values are multiplied by cos^ in order to correct 
the light output for track inclination with respect to direction perpendicular 
to the scintillator segment area. For each scintillator segment in the FRH 
we plot the values of ADC ■ cos 6' as a function of the scattering angle 6 (see 
Fig. 5.1). We parametrize this dependence using a fit with the third order 
polynomial f{0). The ADC values for the analysed events are then corrected 
in the following way: 

ADC- ^ . ,5.1, 




Figure 5.1: 'Raw' ADC value multiplied by cos(^) as a function of 9 for the 
third element of the first layer of FRH. The fitted third order polynomial is 
shown solid blue line. 

For determination of the relation between the corrected ADC values and 
the deposited energy we compare some characteristic points visible in the 
experimental ADC spectra with results of simulations. For this, we create 
two dimensional plots of the corrected ADC values for pairs of overlapping 



segments from consecutive FRH layers. Figure p72\ shows an example of such 
spectrum created for a pair of segments from the FRHl and FRH2 layer. 
The characteristic 'banana' bands corresponding to protons, deuterons, ^He 
and "^He ions are visible. The following characteristic points are well defined 
in this plot: 
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• minimum of ionization for protons (1), 

• punch through point for protons (2), 

• punch through point for deuterons (3), 

• punch through point for '^He (4). 

The punch through points correspond to situation where particles are almost 
stopped in the two scintillator layers. 
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Figure 5.2: Distribution of corrected ADC values from one element of the 
first layer of FRH (y-axis) and adjacent element of the second layer of FRH 
(x-axis). The indicated characteristic points (1)...(4) are explained in the 
text. 

The set of the characteristic points is compared to corresponding values 
obtained from the Monte Carlo simulations, and a linear fit is then used 



as a calibration relation (see Fig. 5.3). The deviations of the characteristic 
points from the linear fit are of about 5 MeV and this value was taken as an 
uncertainty of the energy calibration for the first three layers of the FRH. In 
the case of the 4th layer, the punch through point for the ^He ions couldn't 
be determined due to a very poor statistics. Therefore, the calibration, which 
was based on deuteron and proton bands only, is less accurate than the ones 
for the first three layers. 



5.2.2 Scintillating Electromagnetic Calorimeter 

The main purpose of the Scintillating Electromagnetic Calorimeter is mea- 
surement of energies and emission angles of photons and, therefore, the en- 
ergy calibration is performed to optimize the reconstruction of these particles. 
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Figure 5.3: Deposited energy obtained from MC simulations (y-axis) versus 
corrected ADC value for the third element of the first layer of FRH (x-axis). 
The marked points correspond to: (1) minimum ionizing point for protons, 
(2) punch through point for protons, (3) punch through point for deuterons, 
(4) punch through point for ^He. The fitted linear function represents the 
calibration relation. 

The initial calibration of the SEC crystals was done before the installation of 
the WASA detector at COSY, using cosmic muons and radioactive sources 
[50] [72] . The standard calibration method which corrects possible variations 
of the gain of individual SEC detectors is based on reconstruction of the in- 
variant mass of neutral pions 7r° from two detected gammas originating from 
the 7r° — )■ 77 decay. In the calibration process events with exactly two reg- 
istered neutral particles are selected. The particles are regarded as gammas 
and the invariant mass is calculated according to the following formula [71] : 



where:i?^i andi?-y2 are the measured energies of the gammas based on pre- 
liminary calibration constants, p^i and p^i are their momentum vectors, 61^2 
is their opening angle and ki and ^2 are calibration correction factors for two 
crystals with the highest energy deposited by the pair of the gammas. The 
deviation from the actual invariant mass of tt" is used to optimize the values 
of the calibration factors k in an iterative way according the formula: 




(5.2) 



k 




(5.3) 
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This procedure is applied to all crystals in the calorimeter. It corrects for 
differences in the gain of individual detectors as well as for shower losses and 
border effects. 



5.2.3 Straw Tube Chambers 

The calibration of the Straw Tube Chambers FPC and MDC means deter- 
mination of the relation between the measured drift time of electrons and 
the track distance to anode wire. This relation can be derived from the drift 
velocity and it depends on the gas mixture used in the chambers, the voltage 
on the anode wires and the magnetic field. The calibration is performed using 
the so called uniform irradiation method which is based on an assumption, 
that the number of particles dN passing through a single straw tube within 
a distance interval dr from the anode wire depends only on the width of this 
interval and not on the distance to the wire. This can be written as: 

dN = c- dr, (5.4) 

where c is a constant. Division of both sides of the above equation by the 
drift time interval dt corresponding to the drift distance dr gives: 

dN dr . . 

where v{t) is the drift velocity being a function of the drift time t. Integration 
of this equation over drift time interval from t^m to t gives: 

* dN /■* 

— c/t'= / c-v{t')dt' = c[r{t)-r{t^in)]. (5.6) 

The tmin is chosen as a minimum drift time corresponding to drift distance 
equal to zero {r{tmin) = 0). In this way the searched distance-drift time 
relation can be written as: 

r{t) = '^^ . (5.7) 



The constant c can be calculated with Eq. 5.7 using requirement that the 



maximum drift time tmax corresponds to the radius of the straw tube R 



R 



(5. 



Left panel of Fig. 5.4 presents drift time spectrum used for the calibration. 



The middle panel presents result of integration of the drift time spectrum 
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over the time interval from tmin to t according Eq. |5.6| and the right panel 



shows the resulting distance-drift time relation r(t). As one can see from the 
figure, this relation is close to linear dependence except for the region close 
to the wire and the one close to the straw tube wall. 
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Figure 5.4: Drift time spectrum (left), integrated drift time spectrum af- 
ter subtraction of background (middle) and final drift time-distance relation 
(right). The vertical red lines indicate the range of drift time taken for in- 
tegration and the blue horizontal lines indicate the range of integrals of the 
counts. The figure is adopted from ^S^. 



5.3 Analysis of the ^HepTr events 



5.3.1 Scheme of analysis 

A general scheme of the analysis of dd 
Fig. 



— 7- ■^HepvT events is presented in 



5.5 



In the first step, a preselection of events was performed by applica- 
tion of conservative cuts in order to reduce the data sample and to speed up 
the analysis. The next three steps aimed at selection of the dd — )• ^Hepvr" 
events and included (i) the identification of the ^He, (ii) selection of three 
particles final states and (iii) p, vr" identification. For reconstruction of the 
p and TT" momenta which were not measured in the current experiment we 
used the momentum conservation combined with the information about the 
p and 7r~ scattering angles in the LAB system. Selection of events associated 
with the decay of the "^He — r] bound state was performed by setting an upper 
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limit for the CM momenta of '^He ions. Additionally, appropriate cuts based 
on the simulations, were applied to the p and tt~ kinetic energies and the 
p — TT~ opening angle in the CM frame. After normalization to the integrated 
luminosity, these events were used for determination of the excitation func- 
tion which is a basis of our search for the ^He — t] bound state. For study of 
the background processes leading to the ^Hepvr" final state we determined 
the angular and momentum distributions for the final state particles before 
application the above-mentioned cuts on the ^He, p, tt~ kinematic variables. 

The consecutive steps of the analysis are described in the next subsections. 
The final results of the analysis including the excitation function and the 
angular and momentum distributions are presented in the next chapter. 



5.3.2 Preselection of events 

In order to speed up the analysis, a preselection of the raw data was done us- 
ing the following general conditions which are fulfilled by the ^Hepvr" events: 

1. exactly one track from a charged particle in the Forward Detector, 

2. exactly two tracks from charged particles in the Central Detector, 

3. conservative graphical cut on two-dimensional AE—E plot representing 
the energy deposited in the first layer of the Forward Window Counter 
versus energy deposited in all layers of Range Hodoscope in order to 
reduce background from protons and charged pions (see Fig. 5.6), 

4. conservative cut on the energy deposited in the first layer of the FRH: 
Edep(FRHl)e (0.05 - 0.4) GeV as expected for the ^He ions, 

5. scattering angle for tracks reconstructed in the Forward Straw Tracker 
in the range Opo G (3 — 18)° for eliminating tracks laying outside of 
the Forward Detector acceptance. 



5.3.3 ^He identification 

Identification of the ^He ions was carried out by means of the AE — AE 
techniques. According to the simulations of the reaction dd — )■ i^B.eri)hound 
^Hepvr" , in 98% of all cases the ^He are stopped in the second layer of the 
FRH. Therefore, we identify the ^He ions on the basis of the energy losses 
in the first two layers of the FRH. As on can see in Fig. 5/7, a very clean 
separation of '^He and protons and even of ^He ions is possible. 
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Figure 5.5: Sclieme of analysis of dd — > •^Hepvr" events leading to deter- 
mination of excitation function as indicated with blue arrow as well as of 
the angular and momentum distributions of the ejectiles marked with green 
arrows. Details of the presented scheme are explained in the text. 
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Edep(FRH1-FRH5) [GeV] 



Figure 5.6: Comparison of the Monte Carlo simulation (top) and experimen- 
tal spectrum (bottom) of the energy loss in first layer of Forward Window 
Counter (y-axis) combined with energy deposited in all layers of Forward 
Range Hodoscope (x-axis). The preselection region comprising the ^He ions 
is limited by the red solid line. 

5.3.4 Three-body channels 

The requirement of two reconstructed tracks in the Central Detector and one 
reconstructed track in the Forward Detector selects events from three-body 
channels like ^Hep7r~ but it does not fully exclude events with more than 
three particles in the final state where some of particles are not detected. An 
example can be the dd — )■ ^B.epn~7r^ reaction where the neutral pion in the 
final state remains undetected. 

In order to achieve a more restrictive selection of the three-body final 
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Edep(FRH2) [GeV] 

Figure 5.7: Experimental spectrum of the energy losses in first two layers of 
the Forward Range Hodoscope. The selected area for ^He ions is limited by 
the red line. The empty area below 0.05 GeV in the Edep(FRHl) distribution 
is due to the preselection cut. 



states we use the momentum conservation principle which can be written as: 

Pbeam = PHe + Pi + P2, (5.9) 

where pHe denotes the reconstructed momentum of the '^He detected in the 
Forward Detector, and pi, p2 denote the momenta of two charged particles 
detected in the Central Detector. This equation can be rewritten in a slightly 
changed form which defines a new vector p^ being a difference of the beam 
momentum and the ^He momentum vector: 

Pi + P2 = Pbeam - PHe = Pk- (5.10) 

The vectors Pk,Pi,P2 must lay in one plane and we use this requirement as 
a necessary condition for observation of the three-body final state. For this, 
we define a measure of the coplanarity as an angle between the vector and 
the cross product of the vectors pi and P2: 

Ok,lx2 = APk,Pi X P2)- (5.11) 

In the case of coplanarity Ok,ix2 = SO''. An important advantage of the above 
definition is that in order to calculate the angle 6k,ix2, one only needs to know 
the directions of the pi and p2 vectors and not necessarily the magnitudes. 
This is exactly our case, since due to the lack of the magnetic field in the 
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Central Detector during the present experiment, only the directions of the 
Pi and p2 vectors are measured in the MDC 

The upper panel of Fig. 5.8 shows distribution of the angle 9k,ix2 obtained 
in simulations of the dd — )■ •^Hepvr" reaction. The smearing of the angle 
around the expected value of 90° is due to uncertainties of reconstruction of 
particle tracks. The experimental distribution of 6k,ix2 presented in the lower 



part of Fig. 5.8 has much longer tails which presumably originate from final 
states with more than three particles. 
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Figure 5.8: Distribution of the coplanarity observable 9k^ix2 determined for 
the simulations of the dd ^Hepvr" reaction (upper panel) and for the 
experimental data (lower panel). The cut applied to the data - 9k,ix2 £ 
(90 ± 5°) - is marked by the vertical dotted lines. 

Based on the MC simulations we apply a cut 9kAx2 £ (90 ± 5°) which 
reduces the background coming from reactions with more than three parti- 
cles in the final state. It eliminates also some part of three-particle events 
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with the reconstructed vectors being far from coplanarity due to experimen- 
tal uncertainties. The last case is important in the context of the p -tt~ 
momentum determination method which is based on the assumption that 
the vectors are coplanar as it is discussed in the further part of the present 
chapter. Elimination of events which are far from coplanarity essentially 
improves the resolution of the momentum reconstruction. 



5.3.5 p and tt identification 

The standard method of identification of charged particles in the Central De- 
tector is based on the measurement of energy loss in the Plastic Scintillator 
Barrel and of momentum in the MDC {AE — p method). However, during 
the present experiment the cooling of the solenoid was broken and, there- 
fore, this method cannot be applied. Instead, the energy loss in the Plastic 
Scintillator Barrel was combined with the energy deposited in the Electro- 
magnetic Calorimeter to identify protons and pions [AE — E method). In 
the experimental AE — E spectrum pions and protons are well separated 



which is in line with the simulations (see Fig. 5.9). 



5.3.6 p and tt momentum reconstruction 

As it was already mentioned, the failure of the solenoid cooling excluded 
the possibility of a direct momentum determination for the charged particles 
registered in the Central Detector. However, for the analysed reaction dd — )■ 
^Hepvr" the momentum of the pions and protons can be calculated using 
momentum conservation and the information about their emission angles 
measured in the MDC. 

The reconstruction of the p an 7r~ momenta can be reduced to a simple 
geometrical problem of determination of the pi and p2 magnitudes, assuming 
that we know their directions and that the sum pi + p2 is equal to a known 



vector pk (see Fig. 5.10) 



The magnitudes of pi and p2 vectors are calculated using the relations 
following from the geometrical law of sines: 

Pi = \Pk\ X , \P2\ = \Pk\ X ■ , 5.12 

sm(a; -|- p) sm(a -|- p) 



where the angles a and (3 are defined in Fig. 5.10 

In the calculations we implicitly assume that the vectors pi, p2 and pk 
are co-planar. However, due to experimental uncertainties the directions of 
these vectors reconstructed on the basis of the data from the tracking de- 
tectors might deviate from coplanarity. The applied cut described in the 
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Figure 5.9: Comparison of the Monte Carlo simulation (upper plot) and 
the experimental spectrum (bottom plot) of the energy loss in the Plastic 
Scintillator Barrel (x-axis) combined with the energy deposited in the Elec- 
tromagnetic Calorimeter (y-axis). The green and red curve represents the 
applied graphical cut to separate the protons and pions. 



section 5.3.4 eliminates most of the non-coplanar events and therefore im- 
proves the resolution. 

In order to estimate the precision of the presented method of the pion 
and proton momentum reconstruction, we applied it to simulated dd — )■ 
C^Her]) bound — ^ ^Hepn' events. Fig. 5.11 presents difference between the re- 
constructed and generated momenta. The obtained distribution is symmetric 
around zero, which means that the applied method does not introduce any 
shift in the reconstructed momenta, and the resolution can be estimated from 
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Figure 5.10: Relation between vectors and angles used in the reconstruction 
of the 7r~ and p momenta. 



the RMS value as equal to about 50MeV/c in the LAB frame. 



5.3.7 Cut on ^He momentum 



In the reaction dd — )■ 
spectator 



see 



^He — ri)bound — ^ ^Hepn^ the ^He plays the role of a 
4.1) and is expected to move with low momenta in the CM 
frame corresponding to the Fermi momenta in the ^He nucleus. Hence, we use 
this feature to make a distinction between two regions in the CM momentum 
distribution: the low momentum region (up to 0.3 GeV), where the signal 
from the bound state decay is expected and the high momentum region where 
the direct production dd — )■ ^Hepvr 

The experimental distribution of the 



dominates (see Fig. 5.12). 

^He momenta in the CM system 



is shown in Fig. 5.13 



In the region below 0.3GeV/c, we see a monotonic 
rise of the counts as a function of the momentum. If the process dd — )■ 
C^Her]) bound — B.epn~ was dominating in this region then a structure similar 
to the Fermi momentum distribution should appear with a maximum at 
around 0.2 GeV as visible for the Monte Carlo simulations. 



5.3.8 Cut on p and tt opening angle and kinetic ener- 
gies 



Fig. 5.14 shows experimental distribution of the p — n opening angle in CM 
frame Q^^i before and after the selection of low momentum region of the 
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Figure 5.11: Distribution of reconstructed versus true values of p and vr 
momenta obtained for simulated events (top panel). Difference between re- 
constructed and true values of p and vr momenta in the LAB system (bottom 
panel) . 



■^He ions. The choice of the low momentum ^He is strongly correlated with 
selection of high opening angles 0^_l^^. In spite of it, we apply a cut on the 
opening angle corresponding to the region of (140° — 180°) as expected on 
the basis of the simulations (see Fig. 4.3). This cut removes a small amount 
of events with the opening angle below 140°. 

An additional cut, which was used to improve selection of events from the 
dd — )■ i^B.eri)bound B.epn~ reaction, was applied to the p and vr" kinetic 
energies in the CM system. These kinetic energies are released in the process 
of the 77-meson conversion on a neutron leading to creation of the p — vr^ pair 
and result from the mass difference — m,r- The kinetic energy spectra 
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Figure 5.12: Distribution of the ^He momentum seen in the CM system 
as simulated for the processes leading to the creation of the "^B.er] bound 
state: dd — )■ i^Her]) bound ^B.epn~ (red area) and for the direct production 
dd — Hepvr" (solid line). The blue dashed line shows the selection cut 
applied in our analysis. 




Figure 5.13: Experimental distribution of the '^He momentum seen in the 
CM system. The red dashed line shows the selection cut. 



of protons originating from the bound state decay and from the direct pro- 
duction are similar and are also close to the experimental distributions (see 
Fig. 5.15 -upper plot). Therefore, the applied cut: Ekin^^'^ < 200 MeV is 
not emcient in reduction of experimental background. Larger differences are 
observed for the distribution of the pion kinetic energy which has a maxi- 
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Figure 5.14: Experimental distribution of tlie p — 7r~ opening angle in 
CM frame before (white area) and after (red area) the selection of the low 
momentum ^He. 



mum at around 300 MeV for the bound state simulations compared to about 
220 MeV in the case of the direct production (see Fig. 5. 15 -lower plot). The 
applied cut: Ekin^¥ G (180,400) MeV rejects more than 50% of experimen- 
tal events lying outside the selected region. 



5.4 Luminosity determination 

The present measurements were performed using the ramped beam tech- 
nique. During an acceleration cycle the luminosity could vary due to beam 
losses caused by the interaction with the target and with the rest gas in the 
accelerator beam line, as well as due to displacement of the beam at the target 
position correlated with the variation of the momentum. For this reason it 
was necessary to determine not only the total integrated luminosity but also 
its dependence on the beam momentum. A precise knowledge of this depen- 
dence is important for an accurate relative normalization of data points from 
the studied excitation function of the reaction dd — )■ ^Hepvr" , since too large 
normalization uncertainties can mask the searched signal from the decay of 
the ^He — rj bound state. It is also essential, that the ^He — r] bound state 
does not decay or decays with negligible probability in the channel used for 
the luminosity determination. Otherwise, the normalization procedure can 
cancel the searched signal partly or, in the worst case, fully. 

In the present experiment, the total integrated luminosity was deter- 
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Figure 5.15: Distribution of kinetic energy of p (top plot) and vr^ (bottom 
plot) in the CM frame obtained from experiment and from the MC sim- 
ulations. The red dashed line indicates a boundary of the applied cuts: 
{Ekin^^ < 200 MeV, Ekin^^'^ e (180,400) MeV). 
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mined by comparison of experimental counts from the dd — )■ ^Hen reaction 
witli corresponding cross-sections taken from literature. Because of the lim- 
ited statistics of events registered for this reaction, the relative normalization 
of luminosities corresponding to individual beam momentum intervals, was 
based on the quasi-elastic proton-proton scattering which was observed with 
essentially higher rate. In our analysis we assumed, that the energy depen- 
dence of the cross-sections for the quasi-elastic proton-proton scattering in 
the deuteron-deuteron collisions is proportional to the elastic proton-proton 
cross-sections. 

In the following three subsections we respectively present: basic formulas 
applied for the luminosity calculations, the procedure for the determination 
of the integrated luminosity and the calculations of integrated luminosities 
for individual beam momentum intervals. 



5.4.1 Basic formulas 

In scattering experiments, number of registered events per unit of time for a 
given process can be expressed as: 

dN , ^ , ^ , , 

— {t) = L{t)-a-e, (5.13) 

where a is a cross-section for the process, e is a detection efficiency and 
L{t) is an instantaneous luminosity. In experiments with a fixed target like 
WASA-at-COSY, the luminosity is equal to the product of areal density of 
the target and the number of incident beam particles per unit time. The de- 
tection efficiency e includes geometrical acceptance of the detector, efficiency 
of detection system for registration of the investigated events and efficiency 
of reconstructing events during the analysis. 

The total number of events N registered during the measurement can 



be obtained by integrating both sides of the Eq. (5.13) over the time of 
measurement: 

N = C-a-e, (5.14) 

where C = J^^ L{t)dt is the integrated luminosity. 

In the case of the differential cross- section ^, the Eq. ( |5.14[ ) can be 
modified to: 

dN , , ^ d(7 , , , , , , 

-(x) = £.-W..W. (5.15) 

The integrated luminosity can be determined on the basis of the measured 
counts ^ for a process for which the differential cross-section is known: 
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For the luminosity determination we use angular distributions of the 
cross-section and, therefore, we substitute x = cos 9, dx = d{cos9) = dQ/2iT. 
The previous equation takes the following form: 

£ = ^ i ('5 17) 

^(COS^) ■£(cOSe)' ^ ■ ^ 

The present measurements were performed with a ramped beam and the 
integrated luminosity was calculated for each beam momentum interval sep- 



arately. For this purpose we used Eq. (5.17) reformulated in the following 

way: 



dn 



cos 



Li = / -, (5.18) 

^ (cos 6*, Pi) ■ e{cose, Pi) 

where index i numerates the beam momentum intervals and pi is the central 
momentum for i'th interval. 

For determination of the integrated luminosities £j we use the quasi- 
elastic proton-proton scattering. In our analysis we assume, that the cross- 
section for this process is proportional to the cross-section for the elastic 
proton-proton scattering: 

da da'^^"''^^ 

— {cos e, Pi) = c- ——— {cose, Pi), (5.19) 
d\l d\l 

and c is a proportionality coefficient which can be interpreted in terms of 
shadowing of the protons by the neutrons in the deuteron-deuteron colli- 
sions. The integrated luminosity for the i'th beam momentum interval can 
be written as: 



'^^'(cos^) 1. 



c- ^^^^{cose. Pi) ■ e{cos 6, Pi) c 



The quantities U are proportional to the integrated luminosities and are 
calculated using the elastic proton-proton cross-sections. 

The total integrated luminosity C can be expressed as a sum of partial 
luminosities d over all momentum intervals: 



£ = 5^ A. (5.21) 



The unknown constant c is calculated from the integrated luminosity C which 
was determined on the basis of measurements of the dd — )■ ^Hen reaction: 
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5.4.2 Integrated luminosity 

Cross-section for dd ^Hen reaction 

In order to determine the absolute value of the integrated luminosity we 
used the experimental data on the dd — >■ ^Hen cross-sections measured at 
SATURNE for four beam momenta between 1.65 and 2.49 GeV/c [HI]. As 
shown in Ref. [81] the experimental cross-sections can be very well described 
by a sum of three exponential functions of the four momentum transfer t: 

, 3 

^ = ^a,-e''»(*-*— ), (5.23) 

i=l 

where fit parameters and t^ax denotes the maximum kinemati- 

cal value of t out of the set of measurements for a given beam energy at 
SATURNE. 

In our analysis we used a parametrization of the coefficients Oj, bi with a 
hyperbolical function of the total energy -y/s as proposed in Ref. [82] : 

a^{V^),biiV^) = +r,. (5.24) 

The parameters fcj, g,, obtained from a ffi to the SATURNE data are listed 
in Table 15. 2[ 





ki 


(li 






11.64 


4.05 


-14.49 


bi 


0.78 


3.92 


9.04 


a2 


2327.04 


-1.44 


-3.99.27 


b2 


0.78 


3.92 


9.04 


as 


0.22 


4.08 


1.24 


h 


0.78 


3.92 


9.04 



Table 5.2: Parameters for the total energy ^/s and the four momentum 
transfer t dependence of the cross-section for the reaction dd — )■ ^Hen. In 
the applied parametrization, the a/s and t values are expressed in GeV and 
(GeV/c)^, respectively, and the cross-section is given in /ib/(GeV/c)^. 



The SATURNE data and the applied parametrization are shown in Fig. 5.16 
and [5371 
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Figure 5.16: Total cross-section as function of energy i/s available in the 
CM frame . The solid black line presents the applied parametrization. The 
SATURNE data are marked as black triangles. The red dashed lines indicate 
the range of the energy covered by our experiment {^/s G [4.223, 4.297] GeV). 




Figure 5.17: Differential cross-section as a function of the four-momentum 
transfer. The lines present the applied parametrization. Four sets of the SAT- 
URNE data are plotted as: black full dots (^3=4. 093 GeV), red full squares 
(a/s=4.127 GeV), blue full triangles (a/s=4.159 GeV) and green empty tri- 
angles (v^=4.329 GeV). 
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To determine the integrated luminosity we compared tlie angular distri- 
bution of the measured counts with the parametrized cross-sections in the 
CM system. In this system the four momentum transfer can be expressed 
as: 

t = {Pbeam — PHe) ^ = "^-d^ + mHe^ — 2 • Ebeam " Ehc + 2 ■ \pbeam \ " \pHe \ " COS 9cM , 

(5.25) 

where 6cm is the ^He emission angle in the CM frame. 

To calculate the parametrized cross-section in the CM system we use the 
transformation: 



da 1 da 1 da d{t — t. 



max ) 
~ 1 



dVt 2'Kd{cos 6cm) "^t^ d{t — tmax) d{cos6cM 
where the Jacobian term d{t — tmax) / d{cos 6c m) is equal to: 

d(t tmax 



(5.26) 



d{cos6c 



-'^ ■ \Pbeam\ ■ \PHe\- (5.27) 



CM) 



The statistical uncertainties of the SATURNE cross-sections are much 
smaller than the overall normalization uncertainty equal to 7%. 

Experimental counts 

The measurement of the dd — )■ ^Hen reaction was based on the registration of 
the ^He ions in the Forward Detector. The applied hardware trigger required 
at least one charged particle in the Forward Detector and, in addition, a 
high energy threshold in the Window Counter Detector was set to reduce 
the background from fast protons and pions. At low beam momenta, the 
■^He ions from the dd — )• ^B.en reaction were stopped in the third layer of the 
Forward Range Hodoscope and at high momenta they were stopped in the 
fourth layer. Their identification was based on the AE — E information from 



the Forward Range Hodoscope as it is shown in Fig. |5.18 

The outgoing neutrons were identified using the missing mass technique. 
Fig. 5.19 presents the missing mass of the ^He as a function of the beam 



momentum. Two distinct bands are visible. The lower one corresponds to 
the dd — ^Hen events, whereas the upper one results from pion production in 
processes like dd — ?• ^Henvr^. For the beam momenta above circa 2.3 GeV/c 
the missing mass resolution is worse than at lower momenta. This might 
be due to the fact, that at higher momenta the ^He ions are stopped in the 
fourth layer of the FRH, for which the energy calibration is less accurate 



than for the first three layers (see Chapter 5.2). Also, the visible increase 



of the missing mass as a function of the beam momentum can result from 
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Figure 5.18: Deposited energy in the first layer of tlie Forward Range 
Hodoscope (y-axis) versus deposited energy in all Hodoscope layers (x-axis). 
The black dashed lines indicates the applied cut: Edep{FRHl) > 0.1 GeV 
AND Edep{FRHtot) > 0.3 GeV. 
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Figure 5.19: Missing mass of the ^He as a function of the beam momentum. 



systematic uncertainty of the energy calibration for the fourth layer of the 
FRH. 

For determination of the integrated luminosity we chose the beam momen- 
tum range from 2.239 to 2.293 GeV/c, corresponding to about one quarter 
of the total beam momentum range. On the one hand, the selected interval 
is wide enough to determine the integrated luminosity with statistical accu- 
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racy much smaller than the systematic uncertainty of 7% of the dd — )■ '^Hen 
cross-section (see previous subsection). On the other hand, it is sufficiently 
narrow to neglect uncertainties connected with averaging the cross-section 
as a function of beam momentum, so that for the luminosity determination 



the cross-section for the central momentum can be taken (see Fig. 5.20). 
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Figure 5.20: Parametrization of the dd — )■ ^Hen cross-section as a function 
of cos(^c'Af) for three different beam momenta from the analysed beam mo- 
mentum range: the lowest value pheam = 2.239 GeV/c (red line), the middle 
value 2.266 GeV/c (blue line) and the highest values 2.2293 GeV/c (black 
line). 



Fig. 5.21 shows distribution of the kinetic energy versus the scattering 
angle in the LAB frame for the identified ^He ions. 

The counts for the dd — )• ^Hen reaction cover the angular range from 
about 4° to 10°. As shown in Fig. 5.22 in the CM frame this range corre- 
sponds roughly to the cos 9cm interval from 0.88 to 0.98. For comparison 
of angular distribution of the experimental counts with the cross-sections we 
divided the above cos6'cAf range into five bins of equal width as indicated in 
Fig. [5:221 

The numbers of dd — )■ ^B.en events in the chosen cosOcm intervals were 
obtained from counts in the neutron peak in the corresponding missing mass 
spectra shown in Fig. |5.23[ The number of counts in the neutron peak was 



corrected for the background under the peak approximated with a quadratic 
function which was fitted to a few points selected on the left and on the 
right side of the peak (see Fig. 5.24). The uncertainty of the background 
subtraction was estimated by comparison of fits with a linear and quadratic 
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Figure 5.21: Experimental distribution of kinetic energy versus tlie scatter- 
ing angle in the LAB frame. Tlie simulated kinematic dependence for the 
dd — )■ ^Hera is marked as a colour curve. The red colour indicates the area 
with the cross-section two orders of magnitude higher than in other parts of 
the curve according to the applied parametrization. 




12 



Figure 5.22: Kinematic correspondence between the scattering angle Blab 
and cos^cM for the dd — )■ ^Hen reaction. The dashed lines indicate the 
selected intervals in cos6'ca/- 



function. This uncertainty is equal to 8% and is included as a systematic 
error in the luminosity determination. 

In order to check the applied selection procedure for the dd — )■ '^He?7, 
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Figure 5.23: Missing mass spectra of '^He for five cos 6cm intervals: a) 0.88- 
0.9, b) 0.9-0.92, c) 0.92-0.94, d) 0.94-0.96, e) 0.96-0.98. A clear peak around 
the neutron mass is visible. The second peak around 1.1 GeV/c^ originates 
from pion production reactions. 
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Figure 5.24: Distribution of missing mass of '^He for cosOcm ^ [0.92,0.94]. 
The background was fitted with a linear function (blue dashed line) and a 
quadratic function (red solid line). 
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events, we investigated the response of the Scintillating Electromagnetic 
Calorimeter to the outgoing neutrons. Because of relatively small thickness 
of the SEC from the point of view of the nuclear interactions the probability 
that a neutron will interact in the calorimeter material is only about 40%. In 
order to choose the neutron candidates, in addition to the ^He selection cuts, 
we demanded exactly one neutral particle registered in the Central Detector. 
Besides, we request that the difference between the azimuthal angles for the 
outgoing ^He ion and the neutron is close to 180° as expected for a two body 



reaction: \(pfd — (f)cd — 180° | < 10°. As it is presented in the Fig. 5.25, a clear 
correlation in the scattering angle distribution of the forward and the central 
tracks is visible. The experimental distribution agrees nicely with results of 
simulations. This result confirms correctness of the selection procedure for 
■^He, however due to the poor description of the nuclear processes in the SEC 
the coincidence with the neutrons registered in the SEC was not used for the 
luminosity determination. 



Results for the integrated luminosity 

The efficiency term e{cos6cM) for the analysed reaction was determined on 
the basis Monte Carlo simulations. It turned out, that the calculated effi- 
ciency does not depend on the scattering angle 9cm and is equal to 93%. The 
remaining 7% corresponds to events which couldn't be reconstructed mostly 
due to nuclear interactions of the ^He ions in the FRH. 

The obtained number of counts AA^, as well as the determined luminosi- 
ties for each cos( 



interval, are given in Tab. (5.3) 



cos(6'cm) 


AA^ ± Stat 


Lk[nb ^] ± stat 


0.88-0.9 


37 ±6 


29.5 ±4.9 


0.9-0.92 


43 ±6 


21.6 ±3.3 


0.92-0.94 


93 ±10 


31.7 ±3.3 


0.94-0.96 


143 ± 12 


32.3 ±2.8 


0.96-0.98 


222 ± 15 


34.0 ±2.4 



Table 5.3: Number of counts and integrated luminosity values for different 
cos(^C'm) bins. 



The integrated luminosity was calculated as a weighted average of the 
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Figure 5.25: Scattering angle of charged particle detected in the Forward 
Detector versus scattering angle of neutral particle registered in the Calorime- 
ter. The cuts to select the ^He candidates were applied as explained in the 
text. In addition exactly one neutral particle registered in the Central Detec- 
tor was demanded. Finally the co-planarity condition was applied, namely: 
bcd\ < 10°. The upper plot is the dd — )■ ^Hen MC simulation according 



to the phenomenological parametrization (see 5.4.2). The lower plot shows 
the experimental data. The visible signal from the ^Hen is seen in both plots. 



luminosities determined on the basis of 5 individual cos 6cm intervals: 

=5 



E 



r 1 

1 



(5.28) 



This average integrated luminosity equals 30.72 nb^^. It is marked in the 
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Fig. |5.26|with dashed line. The uncertainty of the integrated luminosity equal 



to 11.54% was calculated by adding in quadrature the statistical uncertainty 
(±4.5%), the systematic error due to the background subtraction (±8%) 
and the normalization uncertainty of the SATURNE data for the dd — )■ ^Hen 
cross-sections |81j (±7%). 
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Figure 5.26: Estimated luminosities Lj for different cos 9cm bins. The sta- 
tistical uncertainties are marked as a vertical bars. The integrated luminosity 
equal to the weighted average is marked as a dashed red line. 



In order to illustrate the capability of the present experiment for deter 
mination of angular distributions of the dd — )• ^B.en cross-section in Fig. 5.27 



we show a comparison between the parametrization of the SATURNE data 
for the beam momentum equal to 2.266 GeV/c and the angular distribu- 
tion based on the present data assuming the integrated luminosity equal to 
the weighted average. Only the statistical uncertainties of the present data 
points are indicated. They can be reduced by extending the measurement 
time. More serious are the systematic errors due to the background subtrac- 
tion. Only careful study of the shape of the background can help to reduce 
this error. 

The above integrated luminosity corresponds to the range of the beam 
momenta from 2.239 to 2.293 GeV/c. In order to calculate the integrated 
luminosity in the full range of the beam momenta which was used in the 
search for the r^-mesic '^He (2.192-2.400 GeV/c) we calculated the correspond- 
ing correction according to the Eq. (5.22) using the relative normalization 
of luminosities determined in various intervals of the beam momentum as 
described in the next section. The total integrated luminosity is equal to 
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Figure 5.27: Angular distribution of tlie dd — ?■ ^Hen cross-section at beam 
momentum of 2.266 GeV/c determined assuming tlie final result for the in- 
tegrated luminosity. The vertical error bars represent the statistical uncer- 
tainties. The cross-section parametrization of the SATURNE data is shown 
with red line. 



117.9 ± 13.6 nb^^. The obtained result is in line with a rough estimation of 
the total integrated luminosity of 170 nb~^, based on the beam intensity and 
the target thickness (see Appendix [B|). 
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5.4.3 Dependence on the beam momentum 
Coincident spectra of two charged particles 

In order to determine the luminosity dependence on the beam momentum 
we used the quasi-elastic proton-proton scattering. In the deuteron-deuteron 
collisions, the protons from the deuteron beam scatter on the protons in the 
deuteron target. The neutrons from the colliding deuterons play the role of 
spectators, which means that before and after the collision they move with 
the Fermi momentum in the CM system of their parent deuteron. We use 
the following notation for this process: dd — )■ ppngpUgp. 

The hardware trigger applied for the registration of the quasi-elastic scat- 
tering required at least one charged particle detected in the Forward Detector 
and another one in the Central Detector. In the analysis, a more selective 
condition of exactly one charged particle in the FD and one particle in the 
CD was set. 

In an ideal case of lack of Fermi motion in the deuteron, the sum of 
momentum vectors of the scattered protons in the LAB system is oriented 
along the beam direction. In this case the azimuthal angles for the protons 
registered in the FD and in the CD fulfil the condition: A0 = (pFD — <Pcd — 
180° = 0°. The Fermi motion leads to smearing of A0 as shown in the 



upper plot in Fig. 5.28 representing results of simulations of the quasi-elastic 



proton-proton scattering. The lower plot in Fig. 5.28 represents the difference 
of the azimuthal angles in FD and in CD determined for the experimental 
data. The central peak with maximum at A0 = 180° lies on a high level of 
background originating from many body processes. 

For further analysis we set the condition |A0| < 20° as expected for the 
p — p quasi-elastic scattering on the basis of the simulations. 

Events selected according the above conditions are presented on a plot 
showing the correlation between the polar angles Opo and Opo oi particles 
registered in the FD and in the CD, respectively. 

In this plot we identified events corresponding to the following quasi-free 
processes: 

• dd-> ppUspUsp, 

• dd dbPtrisp, 

• dd ^ PbdtUsp, 
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Figure 5.28: Distribution of tlie difference of azimuthal angles A0 = (ppo — 
(pCD — 180° for MC simulation of the quasi-elastic dd — )■ ppUspUgp reaction 
(upper plot) and for the experimental events (lower plot). The applied cut 
is marked by red dashed lines. 



dd — )• dn^UspUsp, 



where the subscripts b and d denote the particles from the beam and from 
the target, respectively. 

The identification of these processes was based on a comparison with sim- 
ulations taking into account the Fermi motion of nucleons inside the deuteron. 
The correlations between the polar angles 9fd and 9c d obtained in simula- 



tions of the studied processes shown in Fig. 5.29 coincide with ones for the 



experimental counts presented in Fig. 5.30 



In Fig. 5.30, besides the events from quasi-free processes discussed above, 
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Figure 5.29: Correlations between the polar angles Op^) and 6cd obtained 
in simulations for the different quasi-elastic processes: 1) dd ^ dd, 2) dd 

Pbdtnsp, 3) dd — )■ dbPtUgp. 4) dd — > c/vr+nspn^p, 5) dd — pprispUgp. 



there is a pronounced concentration of events in the area OpD G [3 — 4°] 
and 6cD ^ [40 — 60°]. We presume, that these events originate from the 
quasi-free scattering of the neutron from the deuteron beam on the proton 
inside the target: dd — >• pnpspUsp. In the FD we register the spectator proton 
from the deuteron beam and in the CD the recoil proton from the n — p 
quasi-elastic scattering. For this process we expect no correlation between 
the azimuthal angles of particles registered in the FD and CD, and indeed, 
in the experimental data we see no correlation between these two angles for 
the discussed events. 

In order to suppress the events which contain the charged pions registered 
in the Central Detector, we applied an additional cut on the energy deposited 



in the Plastic Scintillation Barrel as showed in the Fig. 5.31 



As one could expect, after this cut the dd — )■ dn'^UgpUsp events disappear 
from the plot showing the correlation between the polar angles 6pD and 6c d 
(see Fig. [K32| 

From the reactions discussed above, for the luminosity determination we 
decided to use the quasi-elastic p — p scattering. This process is registered 
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Figure 5.30: Correlations between the polar angles 6fd and 6cd from the 
experiment. The indicated areas correspond to the: a)dd — )■ ppUspUgp — )■ 
dn~^nspnsp,h) dd — )■ pnpspHsp ,c) dd — )■ dhPt^sp and dd — )■ ppUgpUsp, d) dd — >■ 
ppUspUsp e) (id PbdtUsp . 



with relatively high statistics. Although we can not separate it from the 
quasi-elastic d — p scattering, we show in the next subsection, that it is 
absolutely dominating in the angular region 6fd £ 15 — 17°. 

We also, investigated the background from many-body reactions. For this 
we chose the \4>fd — 4>cd — 100° | < 20° area (see Fig. 5.33) and we checked 



the distribution of counts in the resulting Ofd-(^cd plot shown in Fig. 5.34 



The counts in the angular range 9fd € 15 — 17°, where the quasi-elastic p — p 
scattering dominates, is about 20 times smaller than the p — p signal. 



Cross sections for the elastic p — p and d — p scattering 

For estimation of the cross-sections for the quasi-free proton-proton scatter- 
ing we used the data for the elastic proton-proton scattering at the proton 
beam momentum equal to one-half of the momentum of the deuteron beam. 
For the analysed beam momentum range of 2.192-2.400 GeV/c, the cor- 
responding proton momenta change between 1096 and 1.200 GeV/c. The 
cross-sections for the elastic proton-proton scattering are taken from the 
SAID partial- wave analysis |H3]- We used the energy dependent "Current 
Solution". We checked that its predictions, describes very well the exper- 
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Figure 5.31: Distribution of ttie energy loss in tlie Plastic Scintillator Barrel. 
The experimental data is shown with black line. The results from MC simu- 
lation of the reaction dd — )■ dTT^n marked as green area, whereas the 
results from the dd — )■ pprispUsp are marked as blue area. The MC counts are 
arbitrarily normalized. The applied cut {Edep{PS) > 0.004 GeV) is shown 
as a red dashed line. Please note that the y axis is in logarithmic scale. 



imental cross-sections measured close to our energy range by the EDDA 
collaboration [M] . 

To estimate the contributions from the quasi-free d — p scattering, which 
in our measurements can not be distinguished from the quasi-free p — p 
scattering, we assumed, that it is described by the elastic d—p cross-sections. 
For this, we made use of the experimental p — d cross-sections of Booth et 
al. [85] and of Boschitz et al. [86], measured with deuteron target and proton 
beam with momentum of 990 MeV/c and 1196.2 GeV/c, respectively. The 
equivalent momenta for a deuteron beam scattered on a proton target are 
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Figure 5.32: Correlations between the polar angles OpD and Oqd from the 
experiment after the application of the cut Edep{PS) > 0.004 GeV. 
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Figure 5.33: Selection of events in the |0fd — <Pcd — 100°| < 20° range 
marked with red dashed lines for estimation of background from many-body 
reactions. 



equal to 1.979 GeV/c and 2.391 GeV/c, respectively. These two momenta are 
very close to the lower and the upper limit of the present beam momentum 
range. In Fig. 5.35 the experimental d — p cross-sections are compared with 
the elastic p—p cross-sections calculated using the SAID parametrization. For 
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Figure 5.34: Correlations between the polar angles 6fd and 6cd from the 
experiment after the choice of the background area. 



the forward scattering angles of about 6fd = 17°, the d — p cross-sections are 
about 20 times smaller than the p—p cross-sections. Therefore, we neglect the 
contributions from the quasi-elastic d — p scattering in the current analysis. 

We also checked that the cross-sections for the elastic d—d scattering [87] , 
which in our experiment can not be distinguished from the p — p quasi-elastic 
scattering, within the investigated angular region are smaller than the elastic 
p — p cross-sections by a factor of approximately 50. Therefore, we can also 
neglected contributions from the d — d scattering in our analysis. 



Integrated luminosities for beam momentum intervals 

In order to calculate the integrated luminosity for each beam momentum 



interval we used Eq. (5.20). The number of the quasi-elastic p — p events 
was corrected for the detection efficiency determined in the simulations and 
equal to 53%. Also, the prescaling factor of the applied experimental trigger 
equal to was taken into account. The integrated luminosities for the 
individual beam momentum intervals were calculated using the p — p elastic 
cross-sections. The sum of these luminosities was then normalized to the 
absolute value of the integrated luminosity determined on the basis of the 



dd — )■ Hen measurements using Eq. (5.22). The resulting value of the cor 
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Figure 5.35: Experimental data for the d—p elastic cross-section for deuteron 
beam momentum of plab = 1979GeV/c (upper plot) and plab = 2391 
GeV/c (lower plot) from [85] and [86], respectively. For comparison, the 
elastic p — p cross-sections calculated using SAID for proton beam momenta 
of 0.989 GeV/c (upper plot) and 1.195 GeV/c (lower plot) are shown with 
solid line. 
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recting factor c equals 0.4 and could be understood as a result of shadowing 
of the proton by the neutron inside the deuteron (shadowing effect) which 
reduces the probability of the quasi-elastic scattering. This result should be 
further investigated in measurements of d — d collisions with the WASA-at- 
COSY making use of the momentum analysis of particles registered in the 
CD for a cleaner identification of the p — p quasi-elastic scattering. 

The integrated luminosities calculated for individual beam momentum 
intervals are shown in Fig. 5.36 (bottom). The statistical uncertainty of each 
point is about 4.7%. 
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Figure 5.36: Differential cross-section for the pp — >■ pp elastic scattering 
at ^ = 17° based on the SAID parametrization (upper plot); the number 
of counts of selected dd pp{nn) events (middle plot) and the integrated 
luminosity calculated for individual beam momentum bins (lower plot). 



Chapter 6 

Discussion of results 



In the present chapter we describe determination of the excitation function 
for the dd — )■ {^'^Grfjijound ~^ ^Hepvr" reaction. A fit of the excitation function 
using the Breit-Wigner curve is used for a search of a signal from decay of the 
^He — 7] bound state. We investigate also the experimental background on 
the basis of angular and momentum distributions of the final state particles. 



6.1 Excitation function 

We constructed two types of excitation function for the dd — )■ ^Hepvr" reac- 
tion. They differ in the selection of the events and in the way of normalizing 
the data points. The first excitation function uses events from the "signal- 
rich" region corresponding to the '^He CM momenta below 0.3GeV/c, as it 



was discussed in Sect. 5.3.7 The counts are plotted as a function of the beam 



momentum as it is shown Fig. 6.1[ top). The obtained function is smooth an 



no clear signal, which could be interpreted as a resonance-like structure, is 
visible. A similar dependence is obtained for events originating from the 
"signal-poor" region corresponding to ^He CM momenta above 0.3GeV/c 



(see Fig. 6.1 middle)). We checked also for possible structures in the differ- 
ence between the discussed functions for the " signal-rich" and " signal-poor" 
region. For this, we multiplied the function for the "signal-poor" region by 
a factor chosen in such a way, that the difference of the two functions for 
the lowest beam momentum bin is equal to zero. This difference is presented 
in Fig. |6.1[ bottom). The obtained dependence is flat and is consistent with 
zero. No resonance structure is visible. 

However, we do not treat this result as a final conclusion of non observa- 
tion of the ^He — r] bound state, since one can apply further cuts to reduce 
the background. Additional cuts on the p and tt" kinetic energy distributions 
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and the p — opening angle in the CM system lead us to the construction 
of a second excitation curve. The CM kinetic energies of protons and pions 
originate from the mass deficit — m,r are around 50 MeV and 350 MeV, 
respectively. We select the kinetic energy of protons smaller than 200 MeV 



and of pions from the interval (180, 400) MeV as it is presented in Sect. 5.3.8 



We apply also a cut on the relative p — tt angle in the CM system in the 



range of (140°-180°) (see 5.3.7). The number of selected events in each beam 
momentum interval is divided by the corresponding integrated luminosity. 
In order to use the Breit-Wigner distribution for the description of a possible 
resonance structure in the excitation function, we translate the beam mo- 
mentum intervals into intervals of the excess energy with respect to the ^He?] 



production threshold. The excitation function is presented in Fig. 6.2 It can 
be well described with a quadratic fit resulting in the chi-squared value per 
degree of freedom of 0.98. In the excitation function we observe no struc- 
ture which could be interpreted as a resonance originating from decay of the 
ry-mesic ^He. 
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HepTT cross-section 



We assume, that a signal from the bound state in the excitation curve de- 
termined as a function of the excess energy Q with respect to the ^He — t] 
threshold, can be described by the Breit-Wigner shape: 

a(g,E^^,r,A) = (6.1) 

{Q - EbeY + (2) 

where Ebe is the binding energy, F is the width and A is the amplitude. The 
value of the Breit-Wigner function for the central energy {Q = Ebe) corre- 
sponds to the maximum cross-section for the decay of the //-mesic ^He into 
the ^Hepvr" channel. As it was shown in the previous section the analysis 
of the excitation functions for the dd — )■ ^Hepvr" reaction does not indicate 
any statistically significant enhancement for energies below the 77 production 
threshold. Therefore, we can only determine an upper limit for the cross- 
section for formation of the ^He— bound state and its decay into the ^Hepvr" 
channel. For this, we fitted the excitation function with quadratic function 
describing the background combined with the Breit-Wigner function. In the 
fit we adjusted the quadratic background and the amplitude A of the Breit- 
Wigner distribution. The binding energy Ebe and the width F were fixed 
during the fit. An example of the fit with i^^e—lOMeV and F=10MeV is 



6.2 Upper limit for the dd — )■ {'^Her]) hound ^HepTr cross-section 95 



C 300 
3 

o 
o 

200 



2.2 



C 
3 

o 

u 



1200 — 
1000 



m 800 

c 

3 I 
3 600 — 



400^ 
200 — 

0^ 



200 — 
150^ 
100^ 
50^ 

-50 — 
-100 H 
-150^ 
-200 



2.2 



2.25 



2.3 



2.35 



2.4 



P.O.. 



2.25 



2.3 




p [GeV/c] 

beam 



Figure 6.1: Excitation function for the dd — )■ ^Hep7r~ reaction for the "signal- 
rich" region corresponding to ^He momentum below 0.3 GeV/c (upper panel) 
and the "signal-poor" region with ^He momentum above 0.3GeV/c (mid- 
dle panel). Difference of the excitation functions for the "signal-rich" and 
" signal-poor" regions after the normalization to the lowest beam momentum 
bin is shown in the lower panel. The black solid line represents a straight line 
fit. The beam momentum corresponding to the ^He — rj threshold is marked 
by the vertical red line. 
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Figure 6.2: Excitation function for tlie dd — )• ^Hepvr^ reaction obtained by 
normalizing the events selected in individual excess energy intervals by the 
corresponding integrated luminosities. The solid line represents a fit with 
second order polynomial combined with a Breit-Wigner function with fixed 
binding energy and width equal to -10 and 10 MeV, respectively. The dotted 
line corresponds to the contribution from the second order polynomial in the 
performed fit. The a values are not corrected for acceptance and efficiency 
cuts. 




shown in Fig. 6.2 The fit was performed for various values of the binding 
energy and the width representing different hypothesis of the bound state 
properties. Obtained results comprising the amplitude A and its uncertainty 
a A for the binding energies of -10, -20, -30 MeV and the widths of 10, 20, 



30 MeV, are gathered in the Table 6.1 In each case, the value of the ampli- 



tude A is consistent with zero within the uncertainty a a, which conffrms the 
null hypothesis of non-observation of the signal. 

In order to calculate the cross-sections corresponding to obtained values 
of the amplitude A we still need to take into account the detector acceptance 
for the process dd — )■ {^Her]) bound — ^ ^Hepvr". We performed simulations 
allowing to determine the acceptance as a function of the beam momentum. 
The study of the beam momentum dependence is of high importance, in 
order to exclude the possibility of creation artiffcial signals in the excitation 
curve due to variation of the acceptance. Obtained geometrical acceptance is 
about 60% and the full efficiency including all cuts applied in the analysis is 
about 19% and it varies by only about 1% along the whole beam momentum 



range as it is shown in Fig. |6.3[ 

In order to calculate an upper limit for the dd — )■ {'^B.erj) bound ^ '^HepTr" 
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Table 6.1: Amplitude of Breit-Wigner function and its uncertainty obtained 
from fit with different fixed values of binding energy Ebe and width F. The 
last column contains the maximum cross section for the dd — )■ {^^Qrj)bound — ^ 



'Hepvr process calculated with Eq. 6.2 



cross-section, the oa values obtained in the above described fit have to be 
corrected for the efficiency e (equal to 19%) and multiplied by the statistical 
factor k equal to 1.64485 corresponding to the probability confidence level of 
90%: 

(^■max = • (6.2) 

The obtained upper limits are given in the last column of Table |6.1[ 
One can notice that these limits depend mainly on the width of the bound 
state and not on the binding energy. Therefore, we take average of dmax 
corresponding to different binding energies and we obtain the upper limits 
for the cross-sections of 28, 32 and 41 nb for production and decays of the 
bound state with a width of 10, 20 and 30 MeV, respectively. 

Predictions of the cross-section for creation of the 77-mesic nucleus is a 
very difficult task since it involves computations of behaviour of the many- 
body system and, besides, it depends on the strength of the relatively poorly 
known rj — N interaction. A theoretical model for formation of ^He — 77 bound 
state in the d + d collisions is at present being developed by Hirenzaki et 
al. [88], however, quantitative predictions for the production cross-sections 
are not available yet. 

In order to estimate if the present upper limit for the dd — )■ {^^erf)bound — ^ 
■^Hepvr" cross-section constitutes an essential constrain for existence of the 
77-mesic ^He we constructed a simple model of formation and decay of such 
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Figure 6.3: Dependence on the beam momentum of acceptance for the dd — )■ 
{'^B.er]) bound ^Hepvr" reaction. The geometrical acceptance is shown with 
blue points whereas the full efficiency including detection and reconstruction 
efficiency is shown with red points. 



state. In our model we introduce a hypothesis that the probability of the 
production of the t] meson in continuum and its absorption on the helium 
nucleus is of the same order. Based on this hypothesis we assume that the 
cross-section for the creation of '^B.e — r] the bound state in the d + d collisions 
in the maximum of the Breit-Wigner distribution is about 15 nb as measured 
for the dd — )■ ^He?7 reaction in the vicinity of the kinematic threshold (see 



Fig. 2.11 ). Further on, we assume, that the probability of decay of the He— 77 



bound state into the ^Hepvr channel is equal to | " ^ = |- The factor | takes 
into account the fact that there are four possible r] absorption channels: 

• rjp pn^ (corresponding to dd — )■ ^Hp-K^) 

• rjp UTC^ (corresponding to dd — ^Unir^) 

• r]n rni^ (corresponding to dd — t- ^Yiemt^) 

• rjn p-K^ (corresponding to dd — )■ ^Yiep-K^) 

In turns, the factor ^ represents our guess of the probability that the three 
observer nucleons (ppn), in the process of the r] absorption on the neutron 
in ^He, form ^He in the final state. Estimation of this effect needs, in our 
understanding of the problem, projection of the ^He wave function on the 
wave function of the ^He — n pair and inclusion of the pion and proton 
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rescattering on the observer nucleons. Since, this kind of estimations isn't 
at present in our reach, we justify our guess using an analogy between the 
decay of the ^He — rj nuclei and the ^He-hypernuclei. For the later case it 
was observed namely that in the vr" decay channel the decay mode ^He — )• 
vr^p^He is dominant |89j . 




The above considerations imply that the cross-section for the 
dd — )■ i^Heri) hound Hepvr" reaction should be in the order of 2nb. This 
number is by an order of magnitude smaller than the upper limits for this 
process determined in the present experiment. Therefore, the present exper- 
imental result does not exclude the existence of r^-mesic "^He. 

6.3 Angular and momentum distributions for 
background 

The background observed in the experiment limits the sensitivity of the 
search for the ?7-mesic ^He. Therefore, understanding of reactions contribut- 
ing to the background can be important for the conducted search. We per- 
formed a study of the background on the basis of the momentum and angular 
distributions of the outgoing protons, pions and ^He ions in the measured 
process dd — )■ ^Hepvr". For this, we compared the experimental distributions 
with results of simulations of the following three processes: 

• direct production corresponding to the uniform distribution over the 
phase space available for the reaction products, 

• formation and decay of the ^He — rj bound state, 

• reaction dd — ?■ '^HeA^* — ?■ ^Hepvr^ proceeding with excitation of an 
intermediate N* resonance which subsequently decays in the pn^ pair. 

The first two processes were already discussed. In the following, we present 
some details of simulations of the third process. In these simulations, the 
mass of the intermediate A^* resonance was selected assuming the Breit- 
Wigner distribution with the resonance energy Eq = 1535 MeV and the width 
r = 150 MeV in accordance with PDG ^20j. For a given available energy in 
the CM frame ^/s calculated from the beam momentum, the selection of the 
A^* mass is limited by two conditions: 





(6.3) 



tun* > ra^,-^ + nip. 



(6.4) 
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Eq. |6.3| with the equahty condition can be considered as a limiting case in 



which the whole available energy ^/s goes to the creation of A^* and the 
■^He — A^* pair rests in the CM frame. In turns, Eq. 6A expresses the fact 
that the mass of the resonance must be large enough for its decay into the pvr" 
pair. After the selection of the resonance mass the sum of kinetic energies of 
N* and of ^He in the CM system can be determined in a following way: 



Ekin 



rriN* 



(6.5) 



The direction of the ^He momentum vector is selected assuming its isotropic 
distribution in the CM frame. Next, the simulation of the N* decay into a 
p—TT~ pair is performed under assumption of its isotropic angular distribution 
in the N* rest frame. 

One should stress, that for study of the background reactions we use 
momentum and angular distributions of the outgoing protons, pions and 
^He ions constructed before application of cuts on the '^He momenta and 
on the p and vr" kinetic energy and opening angle. Application of these 
cuts to predictions of various simulation models enforces similarities of the 
investigated distributions and makes the studies of background inconclusive. 

The experimental and simulated angular distributions in the LAB frame 



of He, protons and pions are presented in Figs. 6.4 6.5 and 6.6, respectively. 



In turns the momentum distributions in the LAB system for He, protons 



and pions are shown in Figs. |6.7[ |6.8| and |6.9[ respectively. Study of these 
distributions led us to the following conclusions: 



• Simulations of the direct production and of the excitation of the inter- 
mediate N* resonance provide very similar results. One can understand 
it when taking into account that the width of A* is larger than the range 
of excess energy available for the p — tt~ pair in the CM system. 

• Results of simulations of the bound state differ essentially from ones 
for the direct production in the case of the ^He angular and momentum 
distribution and the pion momentum distribution. 

• The experimental spectra are reasonably well described by the direct 
production. A very good agreement is observed in the case of the pion 
momentum distributions. In turns there is a problem with description 
of the experimental angular distributions of protons and pions for the 
LAB angles below 40°. In the experimental spectra, a rapid drop of 
counts for this angular region is visible. We presume, that this could 
result from an inefficiency of the forward end-cap of the Plastic Scin- 
tillator Barrel, which covers the angles below 40°. 
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Figure 6.4: Distribution of 9 scattering angle of ^He in the LAB frame 
obtained from experiment and from the MC simulations. 




Figure 6.5: Distribution of 9 scattering angle of protons in the LAB frame 
obtained from experiment and from the MC simulations. 
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Figure 6.6: Distribution of 9 scattering angle of vr in the LAB frame 
obtained from experiment and from the MC simulations. 




Figure 6.7: Distribution of momentum of ^He in the LAB frame obtained 
from experiment and from the MC simulations. 
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Figure 6.8: Distribution of momentum of p in the LAB frame obtained from 
experiment and from the MC simulations. 




Figure 6.9: Distribution of momentum of vr in the LAB frame obtained 
from experiment and from the MC simulations. 



Chapter 7 

Conclusions and outlook 



This thesis was devoted to a search for the "^He — r] bound state via exclusive 
measurement of the excitation function for the dd — )■ ^Hep7r~ reaction. Theo- 
retical and experimental background of the search for the 77-mesic nuclei, the 
applied method of measurement, the analysis of data and obtained results 
have been presented. 

The measurement was performed with the internal deuteron beam of the 
COSY accelerator scattered on a deuteron target of the pellet type and with 
the WASA-at-COSY detection system applied for registration of the reaction 
products. During the experimental run the momentum of the deuteron beam 
was varied continuously within each acceleration cycle from 2.185 GeV/c to 
2.400 GeV/c, crossing the kinematical threshold for the r] production in the 
dd ^He77 reaction at 2.336 GeV/c. This range of the beam momenta 
corresponds to an interval of the excess energy in the ^He — rj system from 
-51.4 MeV to 22 MeV. 

In the off-line data analysis, events corresponding to decays of the ?7-mesic 
^He were selected using cuts on the ^He momentum, p and tt^ kinetic energies 
and the relative p — 7r~ angle in the center of mass system. The range of 
the applied cuts was inferred from simulations based on assumption that the 
decay of the ^He — rj bound state proceeds via absorption of the r) meson on 
one of neutrons in the ^He nucleus leading to creation of the p — 7r~ pair, and 
that the outgoing '^He plays a role of spectator and it moves with the Fermi 
momentum in the CM frame. The integrated luminosity in the experiment 
was determined using the dd — > ^Hen reaction and equals 117.9 ±13.6 
The relative normalization of points of the dd — > ^HepTr" excitation function 
was based on the quasi-elastic proton-proton scattering. 

The obtained excitation function is well described with a second order 
polynomial an no resonance-like structure, which could originate from decay 
of the ^He— 77 bound state, is visible. A fit of the excitation function using the 
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Breit-Wigner curve combined with a quadratic term was performed and an 
upper limit for the dd — )■ i^B.eri) hound '^Hep7r~ cross-section was determined 
as equal to 28, 32 and 41 nb for the bound state width of 10, 20 and 30 MeV, 
respectively. 

Sensitivity of the present search for r^-mesic "^He is limited by experimen- 
tal background corresponding to dd — )■ ^B.epn~ reaction proceeding without 
formation of the intermediate ^He — rj bound state. The measured angu- 
lar and momentum distributions of the reaction products are close to those 
simulated under the assumption of uniform phase-space distribution. 

We constructed a simple model of formation and decay of the ^He — 1] 
bound state based on assumption that the formation cross section is equal 
to 15 nb as measured for the dd — > ^Herj reaction close to kinematical thresh- 
old. Estimation results in the dd — )• {^iieri)bound — > ^Hepvr^ cross section of 
2nb being by an order of magnitude smaller than the determined experi- 
mental upper limits for this cross section. Therefore, one can conclude that 
the present measurement does not exclude existence of ?7-mesic ^He. Above 
estimation might be simplistic and more realistic theoretical calculations are 
needed. Indeed, the group of Hirenzaki are developing a theoretical model 
for formation of ^He — f] bound state [HE] and quantitative predictions for the 
dd — )■ {'^B.er]) bound — ^ ^Hepvr" are expected to be available soon. 

Also experimental search for the 77-mesic ^He is being continued with 
WASA-at-COSY PUIIHI]- In November 2010 a new two-weeks measurement 
of the dd — )■ ^Hepvr" channel and also of the dd — )■ ^Henvr^ reaction was 
performed. Analysis of the collected data is in progress. It is planned to 
extend the present study to the '^He — 7] system. The search for 77-mesic 
nuclei is part of the research program at J- PARC [35] and ?7'-mesic nuclei 
have being proposed recently to be included in the research program of GSI 



Appendix A 
Pseudoscalar mesons 



In the quark model, mesons are built of quark-antiquark pairs Quarks 
and antiquarks are fermions with spin s = |. Consequently, the total intrinsic 
spin S of the qq state can be either or 1. By convention, quarks have 
positive intristic parity P and antiquarks have negative intristic parity. If L 
is the orbital angular momentum of the qq system, then the parity P of the 
state equals (—1) * (—1)''" = ("1)^'''"'^ where the first term (—1) is related to 
opposite qq intristic parities and (—1)^ term arises from the space invertion 
transformation. The total spin J of meson is the sum of the intrinsic spin 
S and the angular momentum L. Each quark is assigned additive baryon 
number i? = | and antiquarks are assigned baryon number with opposite 
sign i? = — |. This implies that mesons have baryon number equal to 0. By 
convention, the sign of the flavour quantum number assigned to quark is the 
same as its charge Q. It implies that the flavour carried by charged mesons 
has the same sign as its charge e.g. the strangness of is = 1 and the 
third component of the isospin of vr" is J3=— 1. 

Lets consider only the set of the three lightest quarks {u,d,s) and an- 
tiquarks {u,d,s). Combining these fundamental triplets we can construct 
nine possible combinations of qq states which are decomposed into SU(3) 
flavour-octet states and SU(3) fiavour-singlet state. This decomposition is 
symbolically denoted as: 

3®3 = 8©1 (A.l) 

If we consider only the ground state (relative angular momentum L = 0) 
then the total meson spin J is equal to the intristic spin S and the parity P 
equals -1. Based on the values we can distinguish two nonets of SU(3) 

^We used the ' notation to underline the fact that the flavour of quark q and antiquark 
q' could be different. In the further part of this appendix we will omit ' symbol. 
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mesons: pseudoscalar mesons with = O and vector mesons with = 
The quark composition of the nonet of the hght pseudoscalar mesons is 



presented in the Table A.l 
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Table A.l: Quark composition 



of the nonet of the li; 



ght pseudoscalar mesons. 



Appendix B 



Estimate of the luminosity 
from the target density and the 
beam intensity 

Presented in this appendix calculations of the luminosity are based on several 
assumptions concerning the pellet target parameters and the beam proper- 
ties, which can vary significantly during the acceleration cycle and also from 
cycle to cycle. Hence, the result should be treated as a rough estimate with 
an order-of magnitude precision. 

In experiments with a fixed target like the WASA-at-COSY, luminosity 
is equal to the product of areal target density pt and of number of incident 
beam particles per unit time jb. 

L^Pfjb. (B.l) 
The areal density can be expressed in the following form: 

= (B.2) 

where A^^^ is the number of atoms in one pellet and A is an area of lateral 
beam extension at the target position. The value Nt can be calculated from 
the following equation: 

4 n 
Nt^--7r-Rt'-pD-j^-NA, (B.3) 

where Rt is average radius of pellets, p_D is density of pellet material, n 
is number of atoms per one molecule, M is a moll mass and A^^ is the 
Avogadro constant. For the WASA-at-COSY pellets produced of deuterium 
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(D2 molecules), the following values for the areal density calculations are 
assumed: 

Rt = 15 mm, 
pz. = 0.162^, 
n = 2, 

mol 

The number of incident beam particles per unit time ji, can be estimated 
from the relation: 

jb^f- (B.4) 

where / is the revolution frequency of the accelerator and Nb is the number 
of beam particles stored in the accelerator ring. The values of / and iV^ 
were monitored during the experiment, and are about 1.2 MHz and 3.7-10^ 
particles, respectively. 

Taking all the aforementioned relations and numbers into account, we 
obtain a value for the luminosity of about 10^^ cm~^s~^. However, in the 
calculation we have assumed an ideal 100% overlap of beam and pellet. The 
typical overlap factor in the WASA-at-COSY for a pellet rate of 8 kHz is 
about 0.5. In addition, one needs to take into account the accelerator beam 
losses resulting in another factor of 0.5. Hence, the final estimate of the 
luminosity is L ~ 2.7-10^'' cm~^s~^. Taking into account the effective time of 
measurement T = 62760 s, we obtained a value for the integrated luminosity 
C of about 170 nb~^. 
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